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Abstract 

Conformally-invariant curves that appear at critical points in two- 
dimensional statistical mechanics systems, and their fractal geometry 
have received a lot of attention in recent years. On the one hand, 
Schramm W has invented a new rigorous as well as practical calcu- 
lational approach to critical curves, based on a beautiful unification 
of conformal maps and stochastic processes, and by now known as 
Schramm-Loewner evolution (SLE). On the other hand, Duplantier 
[2 13 has applied boundary quantum gravity methods to calculate ex- 
act multifractal exponents associated with critical curves. 

In the first part of this paper I provide a pedagogical introduction 
to SLE. I present mathematical facts from the theory of conformal 
maps and stochastic processes related to SLE. Then I review basic 
properties of SLE and provide practical derivation of various interesting 
quantities related to critical curves, including fractal dimensions and 
crossing probabilities. 

The second part of the paper is devoted to a way of describing 
critical curves using boundary conformal field theory (CFT) in the so- 
called Coulomb gas formalism. This description provides an alternative 
(to quantum gravity) way of obtaining the multifractal spectrum of 
critical curves using only traditional methods of CFT based on free 
bosonic fields. 
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1 Introduction 



The area of two-dimensional (2D) critical phenomena has enjoyed a recent 
breakthrough. A radically new development, referred to as the Schramm (or 
stochastic) Loewner evolution (SLE) has given new tools to study criti- 
cality in 2D, and also provided us with a new interpretation of the traditional 
conformal field theory (CFT) and Coulomb gas approaches. Examples of 
systems described by SLE include familiar statistical models — Ising, Potts, 
0(n) model, polymers, — as well as "geometric" critical phenomena like per- 
colation, self-avoiding random walks, spanning trees and others. The new 
description focuses directly on non-local structures that characterize a given 
system, be it a boundary of an Ising or percolation cluster, or loops in the 
0(n) model. This description uses the fact that all these non-local objects 
become random curves at a critical point, and may be precisely characterized 
by stochastic dynamics of certain conformal maps. 

The SLE approach is complementary to that of CFT, and the new de- 
scription has not only reproduced many of the known results from previous 
approaches, but also gave new results, either conjectured before or unknown 
altogether. It appears that questions that are difficult to pose and/or answer 
within CFT are easy and natural in the SLE framework, and vice versa. 

The SLE approach is very intuitive and transparent using traditional pa- 
radigms of stochastic processes — Brownian motion, diffusion, and the like. 
In spite of all this, SLE is not yet widely known in the physics community 
and deserves more attention and study. One goal of this paper is to give a 
brief introduction to this burgeoning field. 

Another important recent advance (actually preceding the invention of 
SLE) in the study of critical 2D systems, has been the calculation of exact 
multifractal spectrum of critical clusters by Duplantier IHEI, who has in- 
geniously applied methods of boundary quantum gravity (the KPZ formula 
of Ref. |1]). The second goal of this paper is to review Duplantier 's re- 
sults and rederive them using traditional methods of CFT. To this end I to 
connect SLE with CFT in the so called Coulomb gas formulation. In this 
formulation the curves produced by SLE can be viewed as level lines of a 
height function (bosonic field) that fluctuates and is described by a simple 
Gaussain action with some extra terms. This formalism allows to perform a 
very transparent translation between SLE and CFT, and between geometric 
object (curves) and operators and states in the CFT. 

Many reviews of SLE and its applications in physics already exist. They 
are listed in the references section in the end. Refs. [S|-jlU| are geared 
for physicists, and Refs. JJ-JHI for mathematicians. My presentation 
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is for physicists who may want to read original mathematical papers on 
SLE. Therefore, I use mathematical language and notation, explaining and 
illustrating all important terms and ideas with plausible arguments and 
simple calculations. The presentation is not rigorous, but I try to formulate 
all important statements precisely. Another feature of this paper (mainly 
in its second part) is that I assume that readers are familiar with statistical 
mechanics and methods of CFT, some of which are briefly summarized in 
appropriate sections (see Ref. [lOj for a thorough introduction). 

The structure of the paper is as follows. In Section El I first describe 
microscopic origins of critical curves, as they appear in models of statistical 
mechanics defined on 2D lattices. Then I introduce various quantities of 
interest related to critical curves. 

The next Section |21 presents some properties of conformal maps, espe- 
cially those that map the complement of a curve in the upper half plane 
(UHP) to the UHP. These maps can be obtained as solutions of a simple 
differential equation introduced by Loewner. I provide a heuristic derivation 
of this equation and give a few example of explicit solutions. 

In Section I introduce SLE and describe its basic properties (some 
of which are actually derived later) and relation to particular statistical 
mechanics models. 

Section |S1 provides a quick introduction to tools of stochastic analysis, 
the main ones being Ito formula and its consequences. 

Basic properties of SLE, including its phases, locality and restriction, 
are considered and derived using stochastic calculus in Section El 

In Section [3 I give several examples of non-trivial calculations within 
SLE, whose results provide probabilistic and geometric properties of critical 
curves. This section closes the part of the paper devoted to SLE. 

The remaining sections, based on Refs. |2U1 121j . develop an alternative 
way of analyzing critical curves based on CFT methods in the Coulomb gas 
formalism. 

The name "Coulomb gas" refers to a group of techniques that have been 
very fruitfully used to obtain exact critical exponents of various lattice mod- 
els of statistical mechanics (see Ref. |22]) for a review). A similar method 
was introduced by Dotsenko and Fateev to reproduce correlations of 
the minimal models of CFT. The basic ingredient of all these methods is 
a bosonic action for a Gaussian free field. In Section |H1 1 show how lattice 
models can be related to the bosonic action, and how critical curves can be 
created by certain vertex operators. 

Section begins with definitions and properties of harmonic measure 
and its multifractal spectrum. The identification of curve-creating opera- 
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Figure 1: Ising clusters at low (left) and critical (right) temperature. 



tors as vertex operators of Coulomb gas is used then to derive Duplantier's 
results for many multifractal exponents characterizing stochastic geometry 
of critical curves. 

The final Section [TUl brief! v lists topics related to SLE and its connection 
with CFT, that had to be omitted. This section may serve as a (necessarily 
incomplete) guide to the SLE literature. 

2 Critical 2D systems and critical curves 

Many very simple lattice models of statistical mechanics exhibit critical phe- 
nomena characteristic of continuous phase transitions. The prototype of all 
such models is the Ising model which describes the behavior of a collection 
of "spin" variables Si located on sites of a lattice labeled by the index i, and 
taking values ±1. In this paper we will only consider two-dimensional mod- 
els. In the Ising model the energy of the system is given hy H = —JJ2 SiSj, 
where the sum is over all nearest neighbor pairs of sites ij. 

At a finite temperature T various possible spin configurations {S} of the 
system have probabilities given by Gibbs distribution e~^^'^ /Z, where the 
partition function is obtained by summing over all possible configurations: 
Z = Yl{s} Qualitative picture of this model is that at low tem- 

peratures the Z2 symmetry between the up {Si = 1) and down {Si = — 1) 
directions of the spins is spontaneously broken, and the majority of spins 
points, say, up. As the temperature is increased, typical configurations in- 
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Figure 2: Domain walls in a finite Ising system. The boundary conditions 
change at points A and B, forcing a domain wall to go between these points. 
Left figure: zero temperature. Right figure: critical temperature. 

volve small domains or connected clusters of down spins in the sea of up 
spins. The typical size of such clusters — the correlation length — increases 
indefinitely, as the temperature approaches a specific critical value Tc. At 
this temperature the clusters of up and down spins of all possible sizes are 
mixed together, and the whole picture is scale- invariant, see Fig. lU^ The 
cluster boundaries or domain walls at the critical point are fractal curves 
that the SLE focuses on. 

To be slightly more precise, let us consider a system in a simply connected 
region D with the boundary dD, with a very fine lattice inside (essentially, 
we want the lattice spacing to be much smaller than the system size and 
the correlation length at a given temperature), see Fig. [2 We can force 
a domain wall to go between two points A and B on the boundary dD. 
To this purpose let me impose the following boundary conditions. On the 
upper portion of the boundary between A and B we force the spins to be 
up, and on the lower portion — to be down. Then at zero temperature 
there will be exactly one straight domain wall between the points A and B. 
As the temperature increases, the domain wall will wander off the straight 
line, and eventually, at the critical temperature will become a complicated 
fractal curve. These curves will differ between the members of the statistical 
thermal ensemble, and will have particular weights or distribution within the 
ensemble. 

Another prototypical example of a model that exhibits critical behavior 

^The applet that produced these pictures is available at 
http://www.ibiblio.org/e-notes/Perc/contents.htm 
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Figure 3: Critical percolation on a triangular lattice with a boundary con- 
dition that forces a domain wall (percolation hull) between the origin and 
infinity. The figure is borrowed from Ref. jl7j . 

is the site percolation. In this model each site on a lattice is independently 
colored grey with probability p or white with probability 1 — p. For each 
lattice there is a critical value p = Pc such that an infinite connected cluster 
of grey sites appears in the system (that is, for p < Pc all the gray clusters 
are finite). For a triangular lattice Pc is known to be exactly 1/2. A good 
graphical representation of the critical site percolation on a triangular lattice 
is obtained if we replace every lattice point by a hexagon whose vertices lie 
on the dual honeycomb lattice. Then we can again force a domain wall into 
the system by making hexagons grey and white on two adjacent portions 
of the boundary. The Fig. 13 (borrowed from Ref. JZj) shows the upper 
half plane tiled with such hexagons. All the hexagons to the left of the 
origin on the horizontal axes (this is the boundary) are colored grey, and all 
the hexagons to the right are colored white. This produces a domain wall 
separating grey and white hexagons the beginning of which is shown in Fig. 
131 Fig. 0] shows a much bigger system with the same boundary conditions. 

Both the Ising and percolation models can be included in a larger class of 
models, loosely called the loop models, since their partition functions can be 
written as sums over loop configurations £, either on the original or some 
related (the dual or the surrounding) lattice. Mappings between specific 
lattice models and loop models are described in detail in many reviews 
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Figure 4: A long percolation hull. The figure is borrowed from Ref. jl7j . 



[221 1241 051 126j ■ Here I only mention the 0(n) model because of the richness 
of its phase diagram. The model is defined in the simplest way on the 
honeycomb lattice directly in terms of closed loops: 



Here x is the variable related to the temperature, L is the total length of all 
loops and N is their number in the configuration C. 

It is known from various approaches that the 0(n) model has a critical 
point at some value Xc{n) for all n in the range — 2 ^ ra ^ 2. At the critical 
point the mean length of a loop diverges, but loops are dilute in the sense 
that the fraction of the vertices visited by the loops is zero. For x > Xc{n) 
the loops are still critical but now visit a finite fraction of the sites. This 
is called the dense phase of the loop model. Finally, at zero temperature 
(x = oo) the loops go through every point on the lattice, and this is called 
the fully packed phase. 

For some values of the parameter n the 0(n) model is related to other 
known statistical mechanics models: n = 2 corresponds to the XY model, 
the limit n = describes self-avoiding walks or polymers, and n = —2 
corresponds to the so-called loop-erased random walk. The dense phase of 
the 0(n) model is also related to the critical point of the (/-states Potts 
model. The Potts critical point exists for all ^ g ^ 4, and at that point 




(1) 
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the boundaries of the so-cahed Fortuin-Kasteleyn clusters that appear in the 
high-temperature expansion of the g-states Potts model, are essentially the 
same as the loops in the dense phase of the 0(n) model with n = ^Jq. 

In all the models mentioned above, one can choose boundary conditions 
so as to introduce an open curve starting at one point on the boundary of 
a domain and ending at another boundary point. The continuum limit of 
these open curves with fixed ends on the boundary is exactly what is being 
studied using SLE. 

One crucial paradigm in the study of critical phenomena is that of con- 
formal invariance [271 128j . For the critical curve described above this means 
that if we map conformally the region D in which our system is defined 
into another region D' , then the statistical weights or the distribution of the 
critical curves will be invariant under such mapping. In other words, any 
two curves that map into each other will have the same weight in the corre- 
sponding thermal ensembles. Then we can study these curves in a standard 
simple region, which we choose here to be the upper half plane H with point 
A at the origin and point B at infinity. 

In this setup we may ask various questions about the critical curves. 
Some of them are geometric. For example, we may want to know the fractal 
dimension of a critical curve. More generally, we can imagine that the cluster 
surrounded by a critical curve is charged, and then the charge distribution 
on the domain boundary will be very uneven or "lumpy" . This lumpiness is 
characterized by what is known as the spectrum of multifractal exponents. 
More precise definition uses the notion of the harmonic measure of the cluster 
boundary, and is explained in Section |^ 

Another class of possible questions is probabilistic. We may ask about 
the probability that the critical curve between the origin and the infinity in H 
passes to the left of a given point. Another question asks for the probability 
of the critical curve to touch the boundary at certain places in certain order. 
This is related to the so-called crossing probability in percolation that is 
defined as the probability for a connected cluster to span the critical system 
between two disjoint segments of the boundary. Sometimes we are interested 
only in the asymptotic behavior of probabilities of such events for long times 
or large spatial distances. These asymptotic probabilities behave in a power 
law fashion with some universal exponents that need to be found. 

SLE provides an easy way of answering the above geometric and proba- 
bilistic questions and computing the corresponding quantities. In the follow- 
ing sections I will introduce the necessary tools from the theory of conformal 
maps and stochastic processes and will describe some calculations with SLE. 
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3 Conformal maps and Loewner equation 



Consider the upper half plane M = {z : luiz > 0}, with a curve 7 starting at 
the origin on the real axis such that 7 G HU {0}. We parametrize the curve 
by a real variable t £ [0, cxd), and denote a point on 7 as j{t) and (closed) 
segments as 7[ti,i2]- 

A segment 7[0, t] is an example of the so-called hull. A hull X C H is a 
bounded subset of EI such that H \ is simply connected and K = K CiM. 
So, a hull is, essentially, a bounded (but not necessarily connected) set 
bordering on the real line R. By Riemann's mapping theorem (see Ref. 
[29|). for each such hull there is a conformal map gx that maps M \ K to 
H. Since conformal automorphisms of IH are Mobius transformations with 
real coefficients, we can make gx unique by fixing three real parameters. A 
conventional "hydro dynamic" normalization is such that 

lim (gxiz) -z)=0, (2) 

or, equivalently, that near z = 00 the map has the form 

00 

9Kiz) = ^ + Y.^- (3) 

n=l ^ 

If the hull K is located a finite distance away from the origin, then gxiz) is 
regular at 2: = 0. In this situation it is more convenient for some purposes 
(see Sections 15.31 16.41 below) to consider the map ^k{z) = gxiz) — gk{0), 
normalized as 

$i^(0) = 0, ^>i^(oo)=oo, $'^(00) = 1, (4) 

Since the function gxiz) takes real values on the boundary of IH \ iC, 
the coefficients a„ E M. The coefficient oi = a{K) is called the half-plane 
capacity (or simply capacity) of the hull K. For any r > the map gxiz) 
satisfies the scaling relation grK{z) = rgxiz/r), which implies the scaling 
for the capacity 

a{rK) = r^a{K), Vr > 0. (5) 

Thus the capacity has the dimension of area. Geometrically, it is bounded 
above by R^, where R is the radius of the smallest semicircle that completely 
encloses the hull K. 

Conformal maps for hulls can be composed as shown in Figure |S1 Note 
that the mapping region monotonically shrinks under such a composition. 
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and the hulls grow. Their capacities satisfy another important property 
(additivity or composition rule): 

a[AlJg-/{B)) = a{A) + a[B), (6) 

which can be easily checked by composing the conformal maps qa and gs 
in the form of Laurent expansions (jSJ and finding the coefficient oi of the 
composed map. 

Consider now the map S'-yfo,*] = Qt for the hull that is a segment of a curve, 
as in the beginning of this section. We can always choose the parametrization 
for the curve in such a way that 

a(7[0,f]) =2t. (7) 

We will call the parameter t "time" , since the evolution of gt in this variable 
will be of importance. Then it can be shown that the map gt satisfies a very 
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Figure 6: Illustration for the derivation of Loewner equation. 



simple differential equation called Loewner equation j3Up: 

2 

dt9t{z) = -, go{z) = z, (8) 

gt[z) — l;t 

where S^t is a real function that is the image of the tip of the cut 7(t) under 
the map gt'- 

6 = 9t{l{t)). 

Let me give an intuitive derivation of this equation. Suppose that we 
already know the map gt and want to find out what happens during the time 
increment between t and t + dt. Using the composition of maps we write 
gt+dt = dgt o gf. This composition is illustrated in Fig. IHl Under the map 
gt the segment j[t,t + dt] is mapped to a (almost) straight short vertical 
segment beginning at point ^t G ^- Using the additivity property, Eq. (jEJ, 
the capacity of this little segment is 

a{gt{-f[t,t + dt])) = 2dt. (9) 

The corresponding conformal map dgt removing the segment is elementary: 

dgtiw) = 6 + ^/W^WTMi. (10) 

^This equation has a fascinating history. It was invented in 1923 by Karl Lowner 
(who later changed his name to Charles Loewner, see more about him at http://www- 
gap.dcs.st-and.ac.uk/~history/Mathematicians/Loewner.htmI) to partially solve a famous 
conjecture from the theory of univalent functions proposed by Bieberbach in 1916. After 
many partial successes, the conjecture was finally proved by de Branges in 1985. The key 
element of the proof was the same Loewner equation! A very readable account of this 
story and the proof is given in Ref. |31|. 
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Composing this with gt and expanding in small dt we get 

9t+dt{z) = dgt{gt{z)) = 6 + {gt{z) - it? + Ut 
2dt 

~ gt{z) + 



9t{z) -6 

This immediately leads to Loewner equation (jSJ in the limit dt ^ 0. 

There are two ways in which one can think about Loewner equation. 
The first one was just presented: given a curve 7 in the upper half plane, 
we can obtain, at least in principle, the real function in the equation 
by constructing the corresponding conformal maps. The second way is the 
opposite: given a real continuous "driving" function we can plug it into 
Loewner equation and solve it forward in time starting with the initial condi- 
tion go{z) = z. It is known that the solution exists, but does not necessarily 
describe a map from EI cut along a segment of a curve. In some cases the 
hull that corresponds to the solution gf contains two-dimensional regions of 
the upper half plane, as one of the examples below shows. 

In general, the hull generated by the solution of Eq. (jSJ is defined as 
follows. For a given point z £ M, the solution of Eq. © is well defined 
as long as gt{z) — 7^ 0. Thus, we define Tz as the first time r such that 
limfyT {gt{z) — ^t) = 0. For some points in EI the time Tz = 00, meaning 
that at these points the Loewner map is defined for all times. The union of 
all the points z for which ^ t is the hull corresponding to the map gt{z): 

Kt = {zeW:Tz^ t}, 

and its complement Ht = {z £ M : Tz > t} = M \ Kt is the domain of gt, 
that is the set of points for which gt{z) is still defined. 

Another useful notion is that of the trace 7 produced by Loewner equa- 
tion. This is defined as the union points 

7(t) = limo,-i(z + 6), (11) 

where the limit is taken within the upper half plane. Note that the trace 
and the hull are not necessarily the same objects, as we will see in a simple 
example below, and especially in the case of SLE^ for some values of the 
parameter k (see Section [6. '2() . The reason for this is that points may enter 
the growing hull in two different ways. Some of them are added to the trace 
itself, but others are swallowed, or enclosed by the trace "inside" the hull, 
see examples below. 

One can exhibit many explicit solutions of the Loewner equation for 
several forms of the driving function ^t, see Ref. |S2|- I will give here two 
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Figure 7: The sequence of maps for the construction of the function gs{z) 
in Eq. H12|) . The figure is borrowed from Ref. |32j . 



of them as ihustrations. If S^t 
simply 



c is a constant, the solution of Eq. 



IS 



gtiz) = c + ^/iz-cy+4t. 

The corresponding hull is the vertical straight segment between c and c + 
2i\/t. In this case the map gt can be found by elementary means. 

Another straightforward but instructive example described in detail in 
Refs. |321 133j deals with a circular arc of radius r growing in the complex z 
plane from the point r on the real axis towards the point — r. The segment 
of this arc spanning the angle s S [0, vr) is mapped to an interval on the 
imaginary axis [0,ii?s], where Rs = tan(s/2), by the Mobius transformation 
zi = (z — r)/(z + r), and then removed by the transformation from the 
previous example: Z2 = \/ Rs + z\. Further Mobius transformations are 
necessary to satisfy the hydrodynamic normalization @. This leads to the 
mapping 

'as + Z2 



9s{z) 



ai \a 



Z2 



+ 2- 2a: 



(12) 



where = 1 + R^ = 1/ cos'^(s/2). The first three conformal maps in this 
sequence for r = 1 are illustrated in Fig. [7| 

Expanding the function gsiz) near z = oo we find the capacity of the arc 
to be 2t = r^(l — a~^). After the reparametrization of the arc and the map 
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Figure 8: The trace and the hull for a touching event. Here = 2y^ 6(1 — t) 
for t £ (0, 1) and zero elsewhere. Left shows the situation just before touch- 
ing {t — > 1~); right shows the situation after {t > 1). The trace is the thick 
dark line. The hull consists of that line plus the grey area. That area is 
added to the hull at t = 1. Note that there is a continuum of points added 
to the hull at the time of touching, but only one of these, 7(1), is on the 
trace and is not swallowed. The figure is borrowed from Ref. [51] . 

Qs in terms of t we get the solution of Loewner equation 



corresponding to the driving function = SV?^^ — 2t — 2r. The branches of 
the square roots in Eq. (|13jl have to be chosen in such a way that 



Note that at time r = the map gt changes discontinuously. At any 
time before that the hull of the map is the segment of the arc. But exactly 
at t = r the whole region D = {\z\ < r,lmz ^ 0} (the upper half of the 
disc of radius r) is mapped to the point — 2r: all the points in this region 
are swallowed! 

Let me define this notion more rigorously. We say that a point z G EI is 
swallowed if z ^ 7[0,oo) but z G Kt for some t. In other words, swallowed 
points do not lie on the trace, but get enclosed by the trace in the "interior" 
portions of the hull. The time when a point gets swallowed is called the 
swallowing time for this point. 



9t{z) 



{z - r)2 + 2zVr'^ - 2t + {z + r)J{z + - AzVr"^ - 2t 




(13) 




for \z\ ^ r, 
for \z\ < r. 
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The time r = is the swallowing time for the whole region D. At 
this time the hull of the evolution Kr is the closed semi-disc D, while the 
trace is still the semi-circular arc. One may continue the evolution with the 
driving function = for t > and the trace will continue to grow 

as a simple curve from the point —r on the real axis, while the hull will be 
Kt = D U 7[r^/2, t], the union of the semi-disc and the portion of the trace 
grown after time r. Similar swallowing of a region is illustrated in Fig. |H1 
for = 2a/6(1 — t) for t G (0, 1) and zero elsewhere. 

4 Schramm-Loewner evolution 

The remarkable discovery of Schramm HI was that one can study Loewner 
equation Q with random driving functions and in this way obtain all possi- 
ble ensembles of curves with conformally invariant probabilities. Motivated 
by the conformal invariance of interfaces in statistical mechanical models, 
Schramm had argued that the driving function has to be a continuous sta- 
tionary stochastic process with independent increments. This argument is 
well explained in the existing reviews, here I simply indicate the basic idea. 
First, if we want to produce a curve without branching or self-intersections, 
we need to have a continuous input ^t. Next, for the curves to possess a 
conformally-invariant distribution, the corresponding maps have to be com- 
posed of statistically independent infinitesimal maps of the form (|inj) . To- 
gether with the reflection symmetry this leads to essentially unique choice of 
= \^Bt, where k > 0, and Bt is the standard Brownian motion started 
at ^(0) =0 (that is, Wt = dBt/dt is the white noise with unit strength: 
{WtWs) = 6{t — s)). The resulting stochastic Loewner equation 



and the sequence of conformal maps that it produces came to be known 
as SLEk, where SLE stands for stochastic Loewner evolution or Schramm- 
Loewner evolution. 

Notice that after assuming the hydrodynamic normalization and the 
parametrization in terms of the capacity 0, k is the only important param- 
eter of SLE. A we will see shortly, it completely determines the properties 
of SLE, its hulls and traces. 

Often a shifted version of gt{z) in introduced: 



dtgt{z) 



2 



9o{z) = z 



(14) 



gt{z) - y/liBt 



wt{z) = gt{z) - 



wo{z) 



= z. 
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0<K«4 4<K<8 K>8 

Figure 9: The phases of SLE. The figure is borrowed from Ref. |33j . 

This function satisfies the simple Langevin-type equation 

2 _ • 

dtwt{z) = — -- - \/KBt, wo{z) = z. (15) 

wt{z) 

It is the simplicity of this equation together with the powerful methods of the 
theory of stochastic processes that makes SLE a very versatile calculational 
tool. I will show how to do computations with it in the following sections. 
But first let me summarize the most important properties of SLE (some of 
them will be derived later in Section ISJ. These properties are quite non- 
trivial. Some of them have been rigorously formulated and established in 
Refs. 



• First of all, for all values of k one can still define the trace of an SLE 
as the union of points 7(t) = lim Wi'^{z) (and this limit exists). 

Moreover, the trace is a continuous curve staring at 7(0) = 0, reaching 
infinity as t ^ cxd and never crossing itself (self-avoiding). 

• For ^ K ^ 4 an SLE trace 7 is a simple curve (does not have double 
points). In this case the SLE hull coincides with the trace: Kt = 7(0, i], 
and no point in H gets swallowed. 

• For 4 < K < 8 an SLE trace has infinite number of double points. 
The trace sort of "touches" itself and the real axis at every scale. 
Every time such touching occurs, a whole finite region of the plane 
gets swallowed. As time goes on, almost all the points in IHI (except 
the points on the trace) get swallowed. 

• For K ^ 8 the trace is a space-filling curve (a random analog of the 
Peano curve). This means that no point gets swallowed, but all the 
points in EI lie on the trace. 
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Figure 10: Definition of SLE in an arbitrary simply-connected domain D. 
• The fractal dimension of the trace is (proven in Refs. |361 137j l 



The different behaviors of traces and hulls of SLE for different values 
of K may be called phases in analogy with statistical mechanics. These are 
schematically shown in Fig. |2 

A heuristic derivation of some of these properties given below in Section 
iniwill not be rigorous but will still require the use of probabilistic techniques. 
Therefore, in Section [S] I briefly summarize relevant results of stochastic cal- 
culus. Then in Section [7| I give examples of practical calculations with SLE, 
deriving a number of non-trivial critical exponents and scaling functions. 

So far we have defined chordal SLE gt{z) and its traces and hulls only 
in the upper half plane H. We can now map any simply-connected domain 
to H by a conformal transformation F[z). We fix this function uniquely 
by requiring that 



where z and C are two distinct points on the boundary of D. Then, by 
definition, the chordal SLE in D from z to C is the family of maps ht{z) = 
F~^{gt{F{z))), with a possible random time change, see Fig^l The trace 
of the new SLE is 7* = -F'Ht): and the hulls are = F-'^{Kt). 

Conformal invariance of SLE then means that, first of all, a trace 7* in 
the domain D locally looks the same as a trace 7 in H. In particular, it 




for K ^ 8, 
for K ^ 8. 



(16) 



F{z) = 0, 



F(C) = 00 



F'iO = 1, 



(17) 
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Figure 11: Central charge as a function of k. Notice that k and k' = 16/k 
correspond to the same c. For example, for both k = 8/3 and «; = 6 the 
central charge is zero. These values correspond to self-avoiding walks and 
boundaries of percolation clusters. 

has the same fractal dimension. Secondly, various random events (crossings, 
swallowings, etc.) that correspond to each other under the map have 
the same probabilities. The requirement of such conformal invariance was 
crucial in the original definition of SLE. 

Most importantly for applications in statistical mechanics, SLE pro- 
duces conformally-invariant self-avoiding random traces that are statistically 
equivalent to critical curves in statistical mechanics models. It is actually 
very difficult to make this statement precise, and a lot of efforts has gone 
and is going into establishing the correspondence between SLE for different 
values of n with critical points of various lattice models. 

Correlation functions of local quantities at these critical points are de- 
scribed by CFTs with central charges c ^ 1. Bauer and Bernard [HHl argued 
that SLE describes critical curves in all these CFTs. The relation between 
the SLE parameter k and the central charge happens to be 

(8-3.)(>.-6)^ (18) 
This function is plotted in Fig. 1111 It possesses a remarkable duality, namely 



Duplantier |^ has argued that this duality has a geometric meaning. 
Namely, in the language of SLE, for k > 4 an SLE^ hull Kt has the boundary 
dKt (also called external perimeter or frontier), which locally looks like an 
SLE^/ simple curve with fractal dimension df{K') = 1 + 2/k. 



where 




K 



(19) 
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Lattice model 






df{n) 




References 


Loop-erased random walk 


2 


-2 


5/4 




I, 38 


Self-avoiding random walk 


8/3 





4/3 




1391 

i ! 


Ising model 












spin cluster boundaries 


3 


1/2 


11/8 






Dimer tilings 


4 


1 


3/2 




p5 , 40, 


Harmonic explorer 


4 


1 


3/2 




M 


Level lines of Gaussian field 


4 


1 


3/2 




m 


Ising model 












FK cluster boundaries 


16/3 


1/2 


5/3 


11/8 


ESI 


Percolation cluster boundaries 


6 





7/4 


4/3 


1, 43, 44 


Uniform spanning trees 


8 


-2 


2 


5/4 


mi 



Table 1: Some lattice models for which a correspondence with SLE has been 
conjectured or rigorously established. The dash in the df{K') column means 
that the hull and the trace are the same. 



Notice that the central charge vanishes for k = 8/3 and k = 6. These 
values are special from SLE point of view, since for these values the SLE 
hulls possess very special properties called locality and restriction, corre- 
spondingly. We will consider these properties in Sections 16.31 16.41 

Many different arguments, including comparison of critical exponents, 
have lead to correspondences between SLE^ and specific lattice models that 
we summarize in the Table ^ All these models can be related either to 
the critical point or the dense phase of the 0(n) model, or to the critical 
point of the g-states Potts model. The Coulomb gas methods map these 
models to a Gaussian bosonic field theory with coupling constant g and 
with background and screening charges. Within this theory one can identify 
the curve creating operators and establish a relation between g and k. As a 
result, we get the following relations between the parameters n, q, g and k: 

4 

n = — 2cos7r5i, 9 = —, 2 ^ k ^ oo, (20) 

K 

4 

Q = 2 + 2 cos 27r5, 9 = -, 4 ^ k ^ 8. (21) 

Mathematically rigorous formulation of conjectures related to general Potts 
and 0(n) models can be found in Refs. |15j|i351. 

The critical points of 0(n) model correspond to the range 2 ^ k ^ 4 
{1 ^ g ^ 2), while the dense phase is described by k ^ 4 (0 ^ 51 ^ 1). Parts 
of these ranges correspond to negative n. The critical point of the Potts 
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model is the same as the dense phase of the 0(n) model with n = yjq only 
for positive n, which explains the restriction on k in Eq. ^2\\ . All these 
relations are reviewed below in Section |HJ 

Results presented in Section |H1 allow us to arrive at detailed geometric 
description of critical curves by calculating the spectrum of multifractal 
exponents of the harmonic measure. These exponents include and generalize 
the fractal dimension df{K) This is done in Section IHl 

5 Basic results from stochastic calculus 

To analyze SLE and apply it to the study of critical curves, we need to use 
stochastic calculus. This section provides a brief summary of the necessary 
techniques. 

Here we only consider one-dimensional stochastic processes. As an exam- 
ple it is useful to keep in mind a simple diffusion of a particle on an interval 
(a, b) on the real line. Everything trivially generalizes to higher dimensions. 

All the material in this section and much more is very nicely presented 
inRefs. gSlliSl- 

5.1 Stochastic differential equations, Ito integrals and mar- 
tingales 

A stochastic differential equation (SDE) is, essentially, a Langevin equation, 
which mathematicians like to write in terms of differentials: 

dxt = u{xt,t) dt + v{xt,t) dBf. (22) 

Here the first term in the right-hand side is called the drift term, and in the 
second term Bt is the standard Brownian motion (BM) started at i^o = 
(that is, Wt = dBf/dt is the white noise with unit strength). The process 
xt describes a random "trajectory" of a Brownian particle. The simplest 
example is 

dxt = y/tidBt. 

This describes a simple diffusion with the diffusion coefficient k. 

BM is a Gaussian process with independent increments: for any set of 
times ^ ti < t2 < ■ ■ ■ < tk the random variables 

Bti ,Bt2 - Bt^,... , Bt^ - Bt^_^ 
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are independent and normally distributed with zero means and variances 
ti,t2 — ti, etc. We note here that from this definition it can be shown that 
all of the following are standard BMs: 



-Bt, 

Bs+t — Bs, 

0'~^Ba2t, 



Vs,t > 0, 
Va > 0, 



tB 



reflection invariance, 
time homogeneity, 

scaling, (23) 
time inversion. 



One can write the solution of Eq. (|22|1 as 

xt = xo+ / u{xs, s)ds + / v{xs, s) dBg. 
Jo Jo 

The last term here is an Ito integral defined as the limit of finite sums 

^v{Xs,,Si){Bs^+, - Bs,). 



Note the important property that the integrand is taken always at the left 
end of the time interval. Then it is always independent from the increment 
of the BM multiplying the integrand. This means that upon averaging over 
Bt all the terms in the above sum vanish. The same is true then for the Ito 
integral: 



E^ 



v{xs,s)dBs 



^Jo 



0. 



(24) 



Here E^'[. . .] stands for the expectation value, or the average over the real- 
izations of Bf, and the superscript x refers to the initial condition xq = x. 
Similarly, we have the following Ito isometry: 



E^ 



v{xs,s)dBs 











= E^ 


/ v'^{xs, s) ds 






Uo ^ 



The previous two equations imply another very important property of 
the Ito integral, namely, that it is a martingale. Martingale is, essentially, 
a stochastic process Mt which satisfies the following properties: 



E[|Mt|] < cx), 

E[Mt I history of M up to s] 



Vt, 

yt > s. 
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The second condition (which is formahzed using the notion of filtration 
of cj-algebras Ait to describe "the history of M") contains a conditional 
expectation value (see Section l5.6|) . and it means that if we know that at 
some time s the process M has the value Mg, then the expectation value 
of this process in the future at any moment t ^ s is going to be the same 
Ms independently of time t. In particular, the unconditional expectation 
value E[Mt] = Mq is simply given by the initial value of the martingale. 
Martingales necessarily satisfy stochastic differential equations without drift 
terms: 



dMt = v{Mt,t)dBt. (25) 



5.2 Ito formula 



Next we need the so called Ito formula, which describes a change of variables 
in stochastic calculus. Let me formulate it for ID processes first. Suppose 
that we have an Ito stochastic process xt which satisfies the SDE 

dxt = u{xt, t) dt + v{xt,t) dBt, 

where Bt is the standard BM (that is, Wt = dBt/dt is the white noise with 
unit strength). Now we take any "reasonable" function f{x,t) (it should be 
twice continuously differentiable in both arguments) and define yt = f{xt,t). 
Then yt is again an Ito process that satisfies the SDE 



dt oxt 2 dxf 



Mnemonically, we need to expand to first order in time, but to second order 
in Xt- The quantity [dxtY is found using the rules 

{dtf = dt dBt = dBt dt = 0, {dBtf = dt. 

Thus {dxtf = v^{xt,t) dt, and the SDE for yt becomes the Ito formula: 

fdf{xt,t) 1 2, .d'^f{xt,t) .\9fixt,t)\ 

+ v(xt,t)^^p^dBt. (26) 

OXt 

Often one encounters "time-homogeneous" processes when both u{xt) 
and v{xt) do not depend explicitly on time. In this case the process xt 
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is called a diffusion with drift u{xt) and diffusion coefficient f^(xj). For a 
function yt = f{xt) of a diffusion the previous formula slightly simplifies: 

The coefficients here do not depend explicitly on t, which means that in this 
case the process yt is also a diffusion. The differential operator A appearing 
in the first term here is called the generator of the diffusion xj: 

2ft ^ 1 2. ^ d^fjx) , . .df{x) 



2 (ix^ dx 

we can integrate Eq. (P7j) : 

/(xt) = /(xo) + / i/(x,) + fv{xs)^^^ dBs. (28) 
Jo Jo "a^s 

5.3 Stopping times and Dynkin formula 

It is often interesting to study functions of stochastic processes at various 
random times. Such a random time r is called a stopping time if at any 
moment t we can decide whether r < t or not, or, in other words, whether 
r has happened before time t. In our basic example we can consider, for 
example, the escape time, that is, the first time r when the diffusing particle 
leaves the interval (a, 6): 



r{a,b) 



ini[t : xt i {aM- (29) 



In this example in each realization of our random process we are able to say 
whether the particle is in the interval at time t or outside it. Similarly, we 
can consider the first hitting time of a closed set ^ C M: 

TA = inf [t : xt ^ A\. 

Now we evaluate Eq. (|28|) at some stopping time r: 



f{Xr) = /(Xo) + fAfiXs) ds + rv{Xs)^^ dB, 

Jo Jo dxs 



It can be shown that an analog of the martingale property (|24j) holds for Ito 
integrals with limits that are stopping times (this is related to the so-called 
strong Markov property of the BM, which states that even for a stopping 
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time r the increment -Bt+r — -Br is a standard BM independent from Bt for 
t G [0,t]). Then, taking the expectation values on both sides of the last 
equation we get the so-called Dynkin formula: 

E"[/(a^r)] = fix) + [ / ds] . (30) 

we assumed here that xq = x. Also, to really prove this formula, one needs 
to assume that E^[r] < oo. 

The Dynkin formula is extremely useful when we need to find various 
escape probabilities. Let me consider one example in detail. Suppose, we 
have a diffusion xt started at xq = x £ (a, b). Then at the exit time r = T(^a,b) 
(see Eq. (EHI)) the particle can only escape the interval (a, 6) either at the 
point a or at the point b. Then we can ask the question: "What is the 
probability that the escape happens through the point a?" Formally, we 
need to find one of the quantities 

Pa = P[Xr = a], Pfe = P[xr = 6], 

where P[X] denotes the probability of the event X. It is obvious that these 
two probabilities add to one: 

Pa + Pb = 1. (31) 

We will find another equation relating pa and pb using the Dynkin formula. 
To do this, we consider the expectation value 

B^[fiXr)]=Pafia) + Pbfib). 

For any function f{xt) the LHS of this equation is given by the Dynkin 
formula. But if we find a function that satisfies the equation 

Afix) = 0, 

then the Ito formula (|27() implies that f{xt) is a martingale (no drift term 
in the equation), and the Eq. (jnU]) simplifies to Ei^[f{xr)] = /(x), and for 
such a function we get 

Paf{a)+PJ{b) = f{x). 
Combining this with Eq. H31() . we finally find 

_ fix) - fib) _ fia) - fix) 

""/(«)-/(&)' - /(a) - /(6) • ^''^ 
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Since ^ is a linear differential operator, the function f{x) can usually be 
found explicitly. Often it is expressed in terms of the hypergeometric func- 
tion. 

Let us note that to use the formulas ()32|) . we need any non- constant 
zero mode of A. There is a continuum of such solutions parametrized by 
two constants of integration {A is a second order differential operator), but 
both the additive and the multiplicative constants cancel when a zero mode 
is substituted into Eq. 



5.4 Backward and forward Kolmogorov equations 

Let us denote 

b{x,t) = E^[f{xt)]. 

Then taking E^' of both sides in Eq. ()28() and differentiating with respect 
to t, we get 

It turns out that the right hand side here can be expressed in terms of b{x, t) 
also. Roughly speaking (this is not very trivial), the expectation value and 
the operator A can be interchanged (after this A acts on the variable x), 
giving the so called backward Kolmogorov equation: 

§ = i6, b{x,0) = f{x). (33) 

Note that this is different from the more familiar Fokker-Planck equation. 
In fact, the Fokker-Plank equation (called the forward Kolmogorov equation 
in mathematics) involves the operator A* that is adjoint to A: 

A* fix) = \^2i^Hx)f{x)) - ^{U{x)f{x)). 

The forward Kolmogorov equation involving A* appears as follows. The 
process xt has the transition measure density pt{y, x), which means that the 
expectation values of functions of Xf can be found like this: 



E"[/(^*)] = Jf{y)ptiy,x)dy. 

This is equivalent to pt{y,x) = E^[(5(xf — y)], which is a familiar definition 
of the probability density for the process xt- The density pt{y,x) is also 
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the kernel or the Green's function of the diffusion xt- It is this function 
that satisfies the forward Kolmogorov equation with respect to the final 
coordinate y: 

d 

-Q-^Pt{y,x) = A*ypt{y,x). 

Because the operator A* has all the derivatives on the left, the total proba- 
bility is conserved: Jptiy,x) dy = E^[l] = 1. 

5.5 Feynman-Kac formula 

A simple generalization of the backward Kolmogorov equation H33|) is the 
so-called Feynman-Kac (FK) formula. It concerns the expectation value 



c{x,t) = E"[exp ( - j'v{xs)ds^f{xt) 



where V{x) is a continuous function such that the integral in the exponent 
converges as i — > oo, f{x) is as before, and xt is a time-homogeneous Ito 
process (a diffusion). The FK formula is the following partial differential 
equation for c(x,t): 

Be 

— = Ac-Vc, c{x,0)=f{x). (34) 

This formula is obtained (schematically) as follows. We define 

Dt{x)= [v{xs)ds, C(xt,0 = e-^'(")/(x4). (35) 



The process C{xt,t) explicitly depends on t through its first factor, and the 
Ito equation for it is obtained from the formula 

dC{xt, t)=[A- V{xt)] C{xt, t)dt + v{xt)^^^^^ dBt. (36) 
Upon averaging the last term vanishes, as usual, and we get 

-=W[{A-V{x,))C{xut)\. 

Similar to the case of the backward Kolmogorov equation, the right hand 
side can be expressed in terms of c(a;,t), which results in Eq. (PHl. 
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There is a variant of the FK formula that we can call a stationary FK 
formula. Namely, we can choose the function f{x) in Eq. (|35() to satisfy the 
stationary version of Eq. ()34() : 

[A-V{x)]f{x) = 0. (37) 

Then the process C{xt,t) defined in Eq. H35|) with f{x) being a solution of 
H37() is a martingale, since the drift term in the Ito formula (jSHJ vanishes! 
The expectation value c(x, t) is then really a function of x only, and is equal 
to /(x) for all times. 

Now if we know that the process xt is transient, that is, lim^^oo xt = oo, 
we normalize f{x) such that /(oo) = 1, and denote D(x) = D^{x) < oo we 
get 

fix) = hm E^[C{xt,t)] = [e-^(^)] . (38) 

It should be clear now that in this situation we can compute the characteris- 
tic function x{k,x) = E^ [e*'^^^^'-*] of the random variable D. In addition, if 
the variable D{x) is known to be non- negative, the same approach gives the 
Laplace transform L{s,x) of its probability distribution function p{D,x): 

/•CO 

L(s,a;) = E^[e-^^(^)] = / e-'^^''^p{D,x)dD. (39) 

Jo 

Notice that all the quantities D(x), xik,x), L(s,x), and p{D,x) implicitly 
depend on x, the initial value of the random process xt- 

5.6 Conditional probabilities and expectation values 

Sometimes in the study of random variables and stochastic processes it is 
interesting or necessary to restrict the statistical ensemble of realizations to 
a sub-ensemble satisfying a certain condition. This condition may depend 
on the outcome of a certain random event. For example, for a diffusion 
on the real line we may consider only trajectories that always stay on the 
positive semi-axis, or the ones that happen to be on the positive semi-axis 
at a certain time. Such a restriction of an ensemble is called conditioning. 

Within a restricted or conditioned ensemble we can ask for probabilities 
of various events or expectation values of random quantities. These are 
called conditional probabilities and expectation values. In words we can say: 
"What is the probability of an event A given that an event B happened?" 



29 



Such probability is denoted by P [yl | i3] . It is well known in probability 
theory that conditional probabilities are easily calculated by the formula 

P[^lBl = ^l^i^, (40) 

where P[^ and B] is the unconditioned probability that the events A and 
B both happen, and P[-B] is the unconditioned probability that the event 
B happens. Notice that the Eq. (|4U() only makes sense if the event B has 
non-zero probability F[B] > 0. 

Similarly, given that an event B occurs, we may want to find the ex- 
pectation value of a random variable X, denoted E[X | S]. Conditional 
expectation values have many known properties, but there is no general 
explicit formula for them similar to Eq. (|4fljl . 



6 Basic properties of SLE 

In this Section, based mainly on Refs. [351 1471 148j . we consider the basic 
properties of SLE. Some of them have already been mentioned, and they 
will be here illustrated by plausible arguments. These arguments already 
require some calculations typical for SLE. The main idea of most simple 
calculations with SLE is to look at various random events in the physical 
plane, then see what happens at the same time in the mathematical plane. 
Then we choose a simple real function of the SLE process and study the 
values this function assumes during the interesting events. 



6.1 Scaling 

The scaling property of the Brownian motion (|23j) immediately implies the 
scaling for the SLE processes gt{z) and the growing SLE hulls Kf. Namely, 
we have the following stochastic equivalence: 

9t{z) = -ga'itiaz), wt{z) = -Wa^tiaz), in law. (41) 

a a 

By this we mean that the random quantities on both sides of these equa- 
tions have the same probability distribution. Eq. (|41|) is easily derived by 
observing that the SDE for the right hand side contains ^Bg2i as the driv- 
ing function. The scaling for the SLE processes H41|) immediately implies a 
similar scaling for the SLE hulls: 

Kt = -K^2t in law. (42) 
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6.2 Phases on SLE 



The phases of SLE were aheady described above in Section 01 Here we 
provide a crude derivation of the phases and phase transitions between them. 

6.2.1 Transition at k = 4. 

First we discuss the transition at k = 4. To this end we will fix a point 
X G M on the real axis in the physical plane and consider the motion of its 
image xt = wt{x) up to the time when it hits (which may never happen): 

2 ^ 

dxt = — dt — ^J hidBt, xq = x. 

Xt 

In the mathematical plane, we fix points a and h on the real axis so that 

0<a<x<5<oo. 

Let r be the exit time from the interval [a, 6]. Being continuous, the process 
Xt can exit [a, 6] either through a (with probability P^), or through h (with 
probability 1 — Pa). As described in Sec. 15. HI the probability Pa can be found 
if we know a non-constant zero mode f{x) of the generator of diffusion xt, 
see Eq. 

Now if we take the limits a ^ 0, 6 ^ oo, it becomes the probability for 
Xt to hit in a finite time, which is the probability for the point x on the 
physical plane to belong to the hull: 

b^^a^Q f{a) - f{b) 
In general, the order of limits matters here. If it is reversed, 

a^ob^^ f{a) - f{h) 

is the probability for xt to come arbitrarily close to 0, that is, for the hull 
to come arbitrarily close to the boundary. Note that in order to determine 
these probabilities, only the behavior of a zero mode of A at zero and at 
infinity is necessary. 

The generator for the process xt is 

^ _ K 2d 
2 dx'^ X dx 
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A non-constant zero mode of this operator is 

( 1 1 — — 
„/ ^ I X K for K 7^ 4, 

= 1 I I ( A 
1^ log |x| tor K = 4. 

Substituting this into Eqs. 1)431 144p we find that the answer is independent 
of X, thus the probabihty for the huh to touch the boundary 

p _ / for K < 4, 
~ \ 1 for K > 4. 

For K = 4 the order of hmits is important and we find that 

P = 

These formulas clearly exhibit a sort of "phase transition" at k = 4. 



for K < 4, 

1 for K > 4. 



6.2.2 Transition at k = 8. 

This transition is more subtle. To study it, we fix two points < a: < y < cxd 
on the real axis in the physical plane, and compare the times Tx and Ty when 
they enter the growing SLE hull. For k > 4 both these times are finite. 

It happens that for k < 8 there is a finite probability that the points x 
and y are swallowed simultaneously: 'P[tx = Ty] > 0. On the other hand, 
for K ^ 8, with probability one, < Ty In this case the points on the real 
axis are added to the trace sequentially. The same is true for points in H. 

To make these statements plausible (without giving a real proof), let us 
consider Xt = Wt{x), yt = Wt{y), and qt = log|^. By continuity it is clear 
that ^ Xi ^ 7/j ^ oo for all times, so ^ gj ^ oo. 

If X joins the hull before y (that is, Tx < Ty), then = oo, and the 
probability PIqt^ = cxd] > (as well as P[qT^ = 0] > 0). If the points x and 
y join the hull simultaneously, then qt stays finite (bounded) for all times up 
to Tx = Ty, and PIqtx = oo] = 0. So we need to consider the motion of qt- 

Using Ito formula for qt we get 



dqt = —dyt - —dxt - ^{dytf + ^{dxtf 
yt Xt 2yf 2xf 



Here we need to substitute, as usual, 

2 2 

dxt = — dt — y/ndBt, dyt = — dt — y/UdBt. 
Xt yt 
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This gives the following SDE: 

Notice that this equation is not of standard Ito type, since the coefficients 

of the right hand side depend separately on xt and y^, but not on qt- This 
is easily remedied by a trick that is called "a random time change" . 
This time change amounts to consider a new time variable 

Jo xj 

which is a monotonous function of t (since we integrate a positive quantity 
In differential form the time change is 

di = dt/xf. 

We also need to consider the stochastic process 

I 

Notice that 

[dBt)^ = {dBtf/xf = dt/xl = di. 

Therefore, the process = Bf is the standard Brownian motion with re- 
spect to the new time t ! 

In terms of the new variables the SDE for qj: takes the standard Ito form: 

dq^ = - 2) (1 - e-^i)di+^{l - e-'i)dBi. 

The generator of diffusion for this process is 

and we need to find a zero mode of this operator to study the probability 
P[% = oo]. 

The equation Af = can be easily solved by rewriting it as (prime 
denotes the derivative with respect to q) 
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When we integrate this equation, we can drop the integration constants, 
which are inessential, as was explained in the end of Section 15.31 

log /' = (- - 2) log ( sinh I) , /' = ( sinh | 

Since we now only consider k > 4, the function f'{q) exponentially decays 
as (7 ^ 00, and we can choose (ignoring a multiplicative constant) 



( sinh -j " ds. 



When g — > 0, this integral converges at the lower limit when k < 8, and in 
this case we get P[(7ti, = 0] > 0. On the other hand, when k ^ 8, the function 
f{q) diverges as q ^ 0, and P[gri = 0] = 0, implying that P[(7r^ = 00] = 1. 

6.3 Locality 

Many properties of SLE can be discovered by studying how SLE gets per- 
turbed by distortions of the boundary of the domain where it evolves. Such 
distortions can be described by conformal maps. This setting is similar to 
the definition of SLE in an arbitrary simply-connected domain D in Section 
m but there are important differences. 

Specifically, let us consider a usual SLE^ evolving in the upper half plane. 
Consider a hull A located a finite distance away from the origin. Then the 
SLE trace may hit the hull and have a non-zero overlap with its interior. 
Note that this would not happen for SLE defined in the domain H \ ^ as in 
Section^ Let the hitting time of the hull A be ta- 

Next we consider the image of the SLE under the map <I>a removing the 
hull A from the upper half plane and normalized as in Eq. (0]): 

$a(0)=0, $a(oo) = oo, ^>'a(oo) = 1, (45) 

The SLE hulls Kt get mapped to hulls Kt = ^A{Kt), and these can be 
removed from H by the family of maps gt- Though this procedure works for 
any value of k, for simplicity we illustrate it in Fig. El for k ^ 4, in which 
case SLE hulls are the traces. The maps gt are normalized as 

gt{z) = z ^ — + o{z ), z ^ 00, 

where at is the capacity of Kt , and evolve according to the Loewner equation 

dMz) = (46) 

9t{z) - 
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Figure 12: Various maps in the definitions of locality and restriction. Similar 
figure illustrating a commutative diagram of conformal maps first appeared 
inRef. gH]. 



We want now to find the properties of the driving function and see if 
there is a time change that would make it a Brownian motion ^/kB^. If this 
is the case, then gt{z) is the standard SLE^, at least for t ^ta- This implies 
that the original SLE gt is the same as the SLE defined in the domain H \ A 
according to Section^ see Fig. EH and Eq. (|T7jl . Loosely speaking, we can 
say that in this situation the SLE gt does not feel the presence of the hull 
A until Kt hits A. This justifies the name locality for this property. As we 
will see next, this happens only for k = 6, which was rigorously established 
in Ref. [33| (Ref. |1H] contains a simpler proof that we follow here), and 
is in perfect agreement with the statement that SLEg describes the scaling 
limit of critical percolation interfaces. Even on the lattice such interfaces 
are determined locally, since every site is black or white independently of 
the others. 

To find the capacity at and the driving function we notice that the 
hull KtU A can be removed by a different sequence of maps. Namely, we 
first remove the hull Kt of the original SLE by gt- Doing this, we deform 
the hull A into a hull At- Next we remove the hull At by the map ^t = ^At 
normalized the same way as see Fig. 1 121 where an infinitesimal step from 
t to t + (it is also shown. Then we have a commutative diagram of maps, 
meaning that 

-gti'^Aiz)) = '^t{gt{z)). (47) 
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In particular, for the image of the tip of the trace ^{t) we have 

it = gti^Aiim = MMt))) = '^tiCt), (48) 

where = y/^Bt is the usual SLE driving function. 

Prom the derivation of Loewner equation in Section |51 we already know 
that the capacity of the infinitesimal vertical segment gt{'y[t,t + dt]) is 2dt, 
see Eq. @. When this segment (together with the hull At) is mapped to H 
by <I>t, it simply gets rescaled by <I>J(^(). (Note that the map $f is regular 
away from the hull At, and in particular, at the point ^f. Also, since the 
map preserves a portion of the real axis near ^t, its derivative ^t{S,t) > 0). 
Therefore, the scaling property of the capacity, Eq. implies 

dtat = 2^'ti^tf. (49) 

Taking time derivative of Eq. H47|) and using Eqs. ()46l I49|) we get 

, ~ = dtM9t{z)) + <^'tigt{z)) I 

gt{^A{z))-it 9t{z)-^t 

Denoting w = gt{z) and using Eq. H48|) we simplify this to 



dt<^>t{w) 



The right hand side of this equation is non-singular in the limit w ^ (^f To 
see this we expand it in powers of (w — ^j) (using Mathematica) : 

dtM^^) = -3f l'(6) + (^^7^ - ^^(6)) ("^ - 6) + o {{w - Ctf) . 

(50) 

The first term of the expansion gives dt^t{w)\^^^^ = — 3<I>"(^j). 
Finally, using Ito formula for = ^tiCt), we have 

dit = dtMw)L=^dt + mt)d^t + h'li^tmtf 

= <^'^(^t)dCt + {^-3)<l>U^t)dt. 

For K = 6 and only for this value the drift term in the above equation van- 
ishes, and d^t becomes ^/KdBf; after the random time change di = ^[{^.t^'dt, 
proving locality for SLEg. 
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6.4 Restriction 



In this Section we consider the same setup as in the previous one, but only 
for K ^ 4. In this case SLE traces are simple curves, and there is a non-zero 
probability P[7 n ^ = 0] that a trace 7(0, 00) does not intersect the hull A. 
This is also the probability that the hitting time ta = 00. The collection 
of these probabilities for all hulls A located a finite distance away from 
the origin completely characterizes the distribution of the curve 7, since it 
determines the likelihoods for the curve to go through various places. As 
we will show in this Section, these probabilities can be found for SLEg/3 in 
term of the map normalized again as in Eq. (|45|1 . 

Given that an SLE^ trace does not intersect A, we can map the whole 
trace together with H \ A to H by the map $a and ask what the probability 
distribution of the image $^(7) is. If this distribution happens to be the 
same SLE^, we say that SLE^; satisfies the restriction property. Following 
Ref. j4Sj . where conformal restriction was introduced and rigorously studied, 
we will show below that this happens only for k = 8/3. 

We start by studying the rescaling factor ^t{£,t) as a function of time. 
The second term in the expansion (|5n)) leads to 

2 mt) 3 

and then Ito formula gives 

dmt) = dtmw)\^=^,dt + a>'/(6)det + l^t'md^t)' 

The drift term in this equation cannot be removed by any choice of k. 
However, if we apply Ito formula again to Alj:'^^ = ^[{^t)^, we get 

dMf) _dmt) ^h-i[dmt)? 



dM{w)\^=^, = ^q;^-y:'{^t) 



hMj') mt) 2 mt)' 



^m,,,^( i^-^)^^^mt)' (^ (51) 



= mt)'^^' ' V 2 mt)' ' V2 3 J mt 

If we now choose k = 8/3 and h = 5/8, the drift term in the last equation 
vanishes, which implies that M^:^^^'^ is a martingale. Then, on the one hand, 
the expectation value of M^^^^^ is given by the value of this process at t = 0: 
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On the other hand, let us consider the expectation value at the stopping time 
r when the trace 7 hits for the first time either the hull A or the semicircular 
arc Cr of radius R centered at the origin and completely enclosing A. The 
following argument is very similar to the application of Dynkin formula to 
diffusions on an interval in Section [5.31 

At time r we have two options: either 7(t) hits the hull A (in which 
case T = ta) or it hits the arc Cji (in which case ta> t). Then 

E[m(5/8)] = P[r = ta] KA^rJ'^'^ + P[rA > r] KiCr)^'^'^ . 

In the first case the point S^ta hits a "side" of the hull A^-^ (see Fig. [T^ where 
the derivative $'^^(^7-^) vanishes. Thus, only the second option contributes 

to the expectation value E[Mr^'"^^]. Then if we take the radius R to be very 
large, the point becomes very far from the hull Ar, and the derivative 
^t(?t") tends to 1. Since the expectation value of a martingale does not 
depend on time, by taking the limit i? ^ 00 we obtain 

P[-f n A = 0] =P[ta = 00] = $'^(0)^/^ (52) 

(To make this argument rigorous we would need to show that M^^^^^ is what 
is called a local martingale.) 

The beautiful Eq. (|52() can now be used to show that SLE8/3 satisfies 
restriction property. To do this, let us consider two different hulls A and 
B (as in Fig. both a finite distance from the origin, and calculate the 
conditional probability that the image $a(7) of an SLE trace does not in- 
tersect the hull B, given that the original trace 7 does not intersect the hull 
A. This is done using Eq. (|40|) for conditional probabilities: 

r . s , . n Ph^{AU^A(B)) =0] 



p[jnA = 0] 

According to Eq. ()52() the denominator is equal to ^'j^{0)^^^. The hull 
A U $^^(-B) that appears in the numerator in the last equation is removed 
by the composition ^b°^a (see Fig. [SJ, so the probability in the numerator 
is equal to {^b o «'a)'(0)^/^ = $'5($a(0))5/^$^(0)5/8 = ^'^{Of/^^'^{Of/^, 
where we used the product rule for the derivative and the normalization 
$a(0) = 0. Combining all this, we get 

P[*.(.) n B = « I . n ^ = 0] = = *.^,„)Vs. 
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Thus, the image under <I>^ of the subset of the SLEg/s traces that avoid 
the huh A has the same distribution as SLEg/3, which, by definition imphes 
restriction property for SLE8/3. 

The notion of restriction can be apphed to probabihty measures on sets 
of random hulls in the upper half plane that are more general than simple 
curves. These hulls K must be connected, unbounded, and such that KCiR = 
and C \ K is connected. As was proven in Ref. IIB], there is a one- 
parameter family of conformally-invariant measures on such hulls that have 
the restriction property. For all these measures the probabilities of avoiding 
a fixed hull A, as before, are given by 

P[K n A = 0] = '^'^{0)'' , (53) 

where the restriction exponents h ^ 5/8 characterizes a particular probabil- 
ity measure. Restriction measures with exponents h > 5/8 are not realized 
on simple curves. An important example is given by the so called Brownian 
excursions (two-dimensional Brownian motions conditioned to always stay 
in the upper half plane) with filled-in holes which satisfy restriction property 
with h = 1. 

For any value of k in the interval [0,4] other than 8/3, SLE^ does not 
have restriction property, since no value of h makes M^'^^ a martingale, see 
Eq. and the relation is not satisfied. However, the amount by 

which SLEk fails to satisfy restriction property can be quantified. Namely, 
let Sf denote the Schwarzian derivative of the function /: 

f"'{z) Sff'iz) 



f'(z) 2 If 



f'iz) 2\f'iz) 
Then 



$;(6)"exp(-|^*5<l>,(6)ds 



is a martingale (this is an easy consequence of Ito formula and Eq. I|51jl ). if 
we choose 

6-K (8-3k)(k-6) 
h = , Ck = . (54) 

Notice the appearance of the central charge of the CFT corresponding 
to SLE^. This is quite natural, since distortions of the boundary of the 
domain in which a CFT is defined cause changes in its partition function, 
that depend on the central charge. The exponent h in Eq. (|54|) also has a 
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meaning in CFT: it is the dimension hi^2 of a primary operator that creates 
a critical curve when inserted on the boundary, see Eq. H96|) in Section 
ISl Restriction property has been used to relate SLE with CFT in Refs. 

HSISniEIlElESI- 

Finally, let me mention that restriction measures with exponents h > 5/8 
can be constructed by adding certain Brownian "bubbles" (subsets of 2D 
Brownian motions that are closed loops) to an SLE curve with k ^ 4, see 
Ref. |48|. 

7 Calculations with SLE 

In this section I will give detailed examples of calculations of various prob- 
abilities and geometric characteristics of critical curve using SLE. 

7.1 Left passage probability 

This section is adapted from Ref. j54j . 

Let us fix a point z = x + iy £M in the upper half of the physical plane. 
We may ask whether the trace 7 passes to the right or to the left of this 
point. Formally, this is defined using the winding numbers, as follows. We 
can close the curve 7[0, i] by drawing the arc with the radius |7(t)| from 
the tip of the trace to the point |7(t)| on the positive real axis, and then 
draw the straight segment in M to 0. Then the trace 7 passes to the left of 
z if the winding number of the closed curve defined above is 1 for all large 
times t. Since 7 is transient a.s., there is some random time r such that the 
winding number is constant for t £ (r, cxd). This constant is either or 1, 
since 7 does not cross itself. The random time for a given z is either 00 
(if K ^ 4), or the swallowing time of z. Then the trace 7 satisfies 

1 r(-) X /I 4 3 

P[7 passes to the left of z] = - + ^ — ^- 2FA-, -; -; j • 

2 V^r(^ - i)y V2 K 2 y^J 

The idea of the proof of this statement, as well as similar statements 
about crossing probability, is to see what happens in the mathematical plane, 
and relate events in the mathematical plane with some real-valued random 
functions of the SLE process. For this purpose we consider the image of z 
under the shifted function wt = wt{z), and define: 

ut = Rewt, vt = lmwt, Qt = —- 

Vt 
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Then it is almost obvious, that, 7 is to the left of z if and only if 
lim^ y'T^ qt = 00, and 7 is to the right of z if and only if lim^ /-^-z Qt = —00. 
A heuristic argument for k ^ 4 is as follows. In this case the point z is not 
swallowed (r^ = 00). If the trace 7 passes to the left of z, then a particle 
which starts an un-biased isotropic two-dimensional diffusion from z will hit 
MU 7(0, 00) either in [0, 00) or from the right side of 7 with probability one. 
By definition (see Section 1^1 for details), this probability is the harmonic 
measure of 7(0, 00) U [0, 00) from z. It is conformally invariant, which means 
that the harmonic measure in the mathematical UHP of [^j,oo) from gtiz) 
tends to 1 as t — > 00. And this, in turn, means that limj^oo Qt = 00. 

In the case 4 < k < 8 the point z is swallowed at some time which is 
finite with probability 1. At this time the curve 7 closes a loop around z. 
Then the issue is whether the loop is clockwise or counter-clockwise. In the 
first case for times t close to the harmonic measure of 7(0, t) UM is mostly 
concentrated on [0,oo) and the right side of the curve 7(0, t), which implies 
that limty-Tz Qt = 00. For a counter-clockwise loop the same reasoning gives 
limtyr, Qt = -00. 

Next we will use the Ito formula to write the stochastic equation for qt, 
and then the Dynkin formula to find the necessary probability P[limt^oo Qt = 
00]. So, first we have equations for ut and vt which are just real and imagi- 
nary parts of Eq. (jlSf) : 

dut = r, r, dt — \fKdBt, dvt = „ r, dt. (55) 

uf + vf uf + 

The two-dimensional variant of the Ito formula now gives (no explicit time 
dependence) 



dqt = — aut Trdvt = Trdt dBt. 

Vt f i uf + vf Vt 

This does not look like a standard Ito equation, so we redefine the time 
variable and the noise. First, the new time t is introduced as 

tit) = / ^, dt = ^. (56 

Jo vf vf 

This is indeed a time change, since t{t) is a monotonously increasing func- 
tion. Also, notice that as t — > r^, vt ^ sufficiently fast so that t{Tz) = 00. 
In the new time the equation for q^ becomes 

dqr=4^di-^dBt. 
* ql + l Vt 
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Next we set dBf = dBt/vt and note that {dBt)"^ = {dBtf/v1 = dt/v1 = 
dt. This means that B^ = Bt is the standard Brownian motion with respect 
to the new time variable t. The equation for qi now has the standard Ito 
form: 

dqi = -^-^dt - ^/^dB^. (57) 
Next we find the diffusion generator for this equation: 

2^^+g2 + i dq ■ 

If we study this diffusion on an interval (a, 6), where a < q < b, then the 
probability that q^ escapes this interval through the right hand point 
rather then through the left is given by the second of the formulas (jSSI) 
where f{q) should satisfy the equation 



f + J^fi^) = 0. (58) 



This equation has a constant solution, but this solution is not what we 
need, obviously. The other solution is found by straightforward separation 
of variables (and some specific choice of the constants of integration): 

f-q Aj. 

f^^->-lj^TW- '''' 

We expand the integrand in powers of r and integrate term by term: 



fiQ) 



/ dr ^^(4//^ 

-'0 ZTn 



2\m fA /^\ {_r,'2\m 



(4/k) 



ml ^ 2m +1 ml 

m=0 m=0 



Here {a)m denotes the Pochhammer symbol {a)m = r(a + m)/T{a). Using 
this notation we can write 



1 T{m + a) r(a) (a)„ 1 



m + a r(m + l + a) r(l + a) (1 + a)m a (1 + a) 
and, in particular, l/(2m + 1) = (l/2)m/(3/2)m- This gives 

/(^)-^2L. (3/2)^ ^-^^^H2'k'2'-'^ 



(60) 
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Using asymptotics of the hyper geometric function we see that the solu- 
tion f{q) has finite hmits 

lim /(.) = ±^^|^. 

This shows that the considered diffusion of is transient, meaning that 
with a finite probabihty hm^^^ = oo. Thus, we can finahy take hmits 
a — > — oo, 6 — > oo and get the resuh 

P [7 passes to the left of z] = P [ lim qt = 00] 

^ f(x/y)-f{-^) ^1 r(i) X p /I 4 3 _x^^ 

/(oo) - /(-cx)) 2 + 0Fr(i-l)y' ^2'k'2' y^)' 

When K = 2, 8/3, 4 and 8, the last formula simplifies to 

xy aie^z 1 x argz 1 
1-1 — — I 1 — - 

~II2 ' n oil' ' n' 

7r|z|^ TT 2 z\z\ TT 2 

respectively. 

The value 1/2 obtained for k = 8 is somewhat misleading. The point 
is that if K ^ 8, then the curve 7 densely fills the upper half plane, as was 
mentioned in Section 16.21 and goes through every point, not to the left or 
right of it. This is reflected in the fact that for n ^ 8 the function f{q) in 
Eq. 1)59(1 diverges as g — > itcxD. This divergence means that to determine the 
fate of the process as t — > 00, we need to start with a finite interval (a, b) 
{a < x/y <b) and take the limits a — > —00 and 6 — > 00 separately. In both 
cases we find that P[limtyT-2 Qt = ±cxd] = 0, meaning that qi always stays 
bounded. See a related discussion in Section [7. 31 

7.2 Cardy's formula for crossing probability 

The problem is first posed in a rectangle ABCD. We need to find the prob- 
ability that there is a percolation cluster connecting the left side AB and 
the right side CD of the rectangle, where we impose the fixed boundary 
condition {p = 1). Note that from the point of view suggested by SLE, we 
need to consider not the cluster itself, but one of its "boundaries" , upper or 
lower. Let us consider the lower boundary, which in the continuous limit is 
described by SLEg. Then we see that if there is a spanning cluster, then the 
boundary will necessarily start at the point B, and will reach the side CD 
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without touching the upper side AD. In the opposite case, when there is no 

spanning cluster, the boundary will touch AD before touching CD. 

In fact, this reformulation of the problem can be generalized to any 
K > 4, and we will assume this has been done. 

Next we conformally map the rectangle to the upper half plane using 
the Schwarz-Christoffel formula. The direct mapping ^{z) (from rectangle 
to the UHP) is given by an elliptic function, and the inverse mapping — by 
an elliptic integral. Postponing the details until the end of this section, let 
us assume for now that the images of the vertices of the rectangle are 

^{A) = a < 0, <P{B) = 0, $(C) = c> 0, $(D) = oo. (61) 

Since the crossing probability is conformally invariant (as a property of 
SLE), we are now interested in the following question. Since k > 4, both 
points a and c will be swallowed at some finite random times Ta and Tc- 
The crossing probability then is P[tc < Ta], that is, the probability that the 
point c is swallowed before the point a. 

As should be obvious by now, we need to study the motion of the images 
of the points a, b under the Loewner map. In this case it is easier to use the 
original map (before the shift), so we define 

at = 9t{a), Ct = gt{c), n = • 

Ct - at 

The variable rt is normalized to lie between and 1, and we are essentially 
interested in the probability P[ct- = ^(r)] = P[rr = 1], where r is the escape 
time from [0, 1] for rt- 

The calculations are straightforward: 

2 / 2 2 \ 

d{C - at) = di -dt, d{ct - at) = [ dt. 

at-t. ^ct-i at-^y 



Then 



d{i-at) i-at 

dn = 7 rod[ct - at) 

ct-at [Ct- at)^ 

1 1 ^ 2it 



1 - nJ {ct - at)'^ ' Ct-at' 
Again, this SDE is not of the Ito type, and we perform a time change: 

di = dt/ {ct — at)"^, dBt = dBt/{ct — at). 
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Then the process ft = satisfies the Ito equation 

drr = 2(— i — \dt + y/ndBr. 

The generator of diffusion for this process is 

J K d^ /I I \ d 
A=-— + 2(- 



2 dr"^ \r 1 — rJ dr^ 
and its zero mode /(r) is found by simple integrations as before: 

fir) ■ 



(s(l-s)) 



4/k ' 



Since At > 4, the last integral converges both at the lower and the upper 
limits, when r — > 1. As in the previous section, using Eq. this integral 
can be expressed in terms of the Gauss hypergeometric function: 



OO OO 

/(r)= ds s-''^Y.^A/K)J- = r'-'I^Y: 

™ — A * ™ — n 



^™ 1-4/.^ (4/^) 



„ m + 1 — - m! 
m=0 m=0 K 



-r 



1-1 „ti (2-4/k)„ m! 

1 — - \K K K 

K 

At the ends of the interval for diffusion of rt this function takes the values 
/(O) = and /(I) = T'^{1 - 4./k)/T{2 - 8/k). Substituting this into Eq. 
(|32j) with a = 0, 6 = 1, we get the final result 

= Yi^^^Y^^'V, (1. 1 - 1; 2 - 1; .) . (62) 

As usual, here r means the initial value of the process r^, that is, r = 
—a/{c—a). For k = 6 this reduces to Cardy's formula for crossing probability 
for percolation 

Now we can discuss how to map a given rectangle to the UHP. Suppose 
the horizonal and vertical sides of the rectangle have lengths L and L' . It is 
obvious that the crossing probability is invariant under rescaling. Then we 
need to find the (unique) number < /c < 1 (the so called elliptic modulus) 
from the equation 

IJ_ _ K'{k) 
T ~ 2K{ky 
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where K{k) is the complete ehiptic integral of the first kind, and K'{k) = 
i^(Vl — (in the following we simplify these to K,K'). Next we rescale 
the rectangle and place its vertices as follows: 

A = -K + iK', B = -K, C = K, D = K + iK' . 

It is easy to see then that the function 

1 + sn(z, k) 
1 — k sn(z, k) 

maps the interior of our rectangle to the UHP, and its vertices to 

1 — k 2k 
HA) = —, HB) = 0, $(C) = ^, c^(D)=oo. 

/ 1 _ ^\ 2 

Comparing this with Eq. (1611) we obtain r = I ) , which we need to 

V 1 + A;/ 

substitute to Eq. to get the crossing probability for the rectangle. 



7.3 Fractal dimensions of SLE curves 

SLE curves are fractal objects. Their fractal dimension can be estimated by 
the box counting dimension. Namely, we can ask how the number of small 
disks of radius e required to cover an SLE^ curve scales with e: 

where df{n) is the box counting fractal dimension. Strictly speaking, this 
definition is applicable only for finite curves, but it can be applied for any 
segment of a chordal SLE curve, since all the segments should be statistically 
similar. 

The fractal dimension df{K) is related to multifractal exponents of the 
harmonic measure, and can be obtained from them, as explained in Section 
In this section we use a probabilistic approach. 

The dimension df{K) can be estimated in the spirit of Monte-Carlo meth- 
ods by throwing disks of radius e randomly onto the domain containing the 
critical curve, and then counting the fraction of the disks which intersect the 
curve. Alternatively, we can look for the probability that an SLE curve in- 
tersects a given disk. It is clear that this probability should scale as e2~'^/('*) 
(2 here is the dimensionality of the physical plane), and it is this scaling 
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that can be relatively easily obtained from SLE, with the result (rigorously 
established in Refs. 1^11133) see also an earlier discussion in Ref. jSS]) 

=min (1 + ^,2). (63) 

To derive this scaling we need to introduce some notation and properties 
of conformal maps. First, let D be a domain in the complex plain, dD its 
boundary, and z a point inside D. Denote by dist(2;, dD) the Euclidean 
distance between z and the domain boundary. 

If the domain D is mapped conformally to a domain I? by a function 
z = f{z), then the distance between close points z and z + dz gets multiplied 
by a rescaling factor: \dz\ = \f'{z)\\dz\. The same is roughly speaking true 
for finite distances. More precisely, if d = dist{z, dD) and d = dist{z,dD), 
then a corollary to the famous Koebe 1/4 theorem states that 

d Ad d , Ad , 

Let us denote these bounds by d x d/\f'{z)\ and say that both quantities 
are comparable. 

Now we apply this to the Loewner map wt{z) to estimate the limit of 
the distance dt{z) = dist(2;, 7(0, t) U R) between a point z and an SLE curve 
in the physical plane, as the time goes up to the swallowing time (which 
may be intinite). We use the same notation as in section [7.11 and write 
wt{z) = wt = ut + ivt for the image of the point z. In the mathematical 
plane the distance from the image to the boundary is simply Im wt{z) = vt- 
If we introduce the process 

A(.)=log.""*^^)" 



Im Wt{z) 



Eq. ^ gives dt{z) x e-^*("^). Let us find the SDE for Dt{z). First, the 
z-derivative of the basic SLE equation H15() gives 

5* log w[{z) 



wf{z) 

The real part of this equation is 



ni I M 2ReK(-z) o '"t-Ut ,^r■^ 
dt log \w[{z)\ = 1 ^/;|/^ = 2 ' I 65 
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Combining this with the Eq. (|55() for vt, we get 



Thus, Dt{z) increases with t, and to estimate d{z) = dist(z, 7(0, 00) U ^ 
e~^(^) we need to look at 



D{z) = hm Dt{z) = / , , / ^., dt. 



"t ^ "-t 

As in Section mi we change time according to and get 



D{x/y)=Al , (66) 







where the process qi = ui/vi satisfies the SDE (|57|) and has the initial value 
x/y. As we discussed in the end of Section ITHl if k ^ 8 the process qi 
stays bounded as t — > 00. Then the integral in Eq. (|66() diverges, and 
D{x/y) = 00. This immediately gives that d{z) = and 

df{K ^ 8) = 2, (67) 

consistent with the curve 7 densely filling the upper half plane. 

Now consider the case ^ K < 8. Since d{z) x e'-^^^/^), the probability 
P[A(z) ^ e] that the SLE curve intersects the disc of radius e centered at the 
point z is comparable to (scales in the same way with e as) the probability 
F[D{x/y) ^ — loge]. The latter probability can be estimated if we find the 
asymptotics of the probability distribution function p(L', x/y) for D{x/y). 

We expect that the scaling of Y*[d{z) ^ e] with e should not depend on 
the actual position of z. In fact, the SLE scaling property 1)41(1 implies that 
d{x+iy) has the same distribution as , and thus we are free to choose 

the point z anywhere. To simplify the formulas below, we will now take the 
point z to be X + i. Then the process qi starts at qo = x and is transient, 
that is, goes to 00 or —00. In both cases the integral ()66() is convergent 
and non-negative, and we can use the stationary FK formulas (|371 138|) from 
Section 1^31 to find p{D,x) though its Laplace transform L{s,x). Namely, 
L(s,x) should satisfy 

K d'^L 4x dL 4s ^ _ q 

2 ~^ x2 + 1 dx (x2 + 1)2 ~ ■ 
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The change of variables y = / + leads to the hyper geometric equation 



dL 2s ^ 

L = 0. 

dy K 



The solution of this equation normalized as L(s, x = oo) = 1 is 



r(i-a+)r(i 



a- 1 / 1 X 



2 



1 2 //I 2x2 2s 
The inverse Laplace transform give the probability density for D: 

fc+ioo 



1 pc+too 

p{D,x) = ^ / e'^^''^L{s,x)ds. 



The integration contour should lie to the right of all the singularities of 
L{s,x) in the s-plane. If we deform the contour by moving it to the left, it 
will encircle the poles of L{s,x), and for large D(x) the leading behavior of 
p{D, x) will be determined by the pole with the largest real part. 

Let us now find the singularities of L{s, x) given by Eq. (|68|) . Since the 
hypergeometric function 2F1 (a, b; c; x) is an entire function of its parameters 
a,b,c, the only singularities of L{s,x) are in the prefactor in Eq. H68() . 
Gamma functions have poles when their arguments are non-positive integers: 
1/2 — a±{s) = —n, n ^ 0, which gives the poles at real positions 

Sn = —i- ^ 2n — —n . 

" 8 2 

The largest pole is at sq = —1 + At/8, which gives for large D 
Finally, we have the estimate 

/oo 
p{D,x)dD oc e^-'^/^ 
- log e 

which gives 

df{K<8) = l+'^. (69) 
Together with (|H7|) this establishes Eq. 
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7.4 Derivative expectation 



The absolute value of the derivative of the SLE map |t(;^(2:)| and its moments 
are useful quantities. As for any conformal map is the local measure 

of rescaling introduced by the map. For a critical curve described by SLE 
the moments E[|u;j(2;)|''] where /i G M, are also related to the spectrum of 
multifractal exponents of the harmonic measure, as explained in Sectional 
In this context the derivative should be estimated at a certain distance from 
an SLE curve. Alternatively, of all the SLE curves we should choose only 
the ones that pass closer than a certain small distance e from a given point 
z. This is an example of conditioning introduced in Section [5.61 

Finding the conditional expectation value E[|wj(2:)|''; ^ e] at point 
z in the bulk is a difficult problem that has not been solved so far (see, 
however, how similar quantities are calculated in Refs. jlH l,S5j). However, 
each SLE curve starts at the origin in the physical plane, and no conditioning 
is required to estimate the derivative of the SLE map at real point x on the 
boundary. 

Indeed, in this case Eq. (|(i5)) simplifies (since vt = 0) and gives in the 
notation of Section 16.2.11 



According to the Feynman-Kac formula (|34|1 the expectation value c{x, t) = 
E[m(x)|^] satisfies 

K, 2 2h 

dtc(x,t) = -d'icix.t) + -dxc(x,t) 7Tc(x,t), c(x,0) = 1. 

Z X x'^ 

From the SLE scaling law (jSJ we know that x and t must appear in 
the combination x/\/t. With some hindsight we denote y = x^ jlKt and 
c(a;,t) = /(y). Then the equation and the boundary value for /(y) become 



In the limit y ^ (long times) we can neglect y in the brackets in front 
of /', and the equation simplifies to an Euler equation with the solution 




2 



,2 ' 




(70) 



y'r(y) + + ^ + 2/)/'(y) - -/(y) = o, lim /(y) = i 



A{h)/2^ where 



A{h) 



K-4 + - 4)2 + 16k/i 

2k 



(71) 



is a solution of the indicial equation that vanishes as /i — > 0. 
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In the context of the problem of multifractal exponents of harmonic 
measure the scahng 




is all we need. But we can also solve the problem completely. Namely, if we 
write f{y) = e~yy^^^^^'^ilj{y), then ■(/'(?/) satisfies 

yi>"{y) + (c - y)ij'{y) - a^^iy) = 0, V'(y ^ oo) ^ e^y-^^^^^/^ 

2 1 2 1 

a = -H \ c=-H hA(/i), 

which is the standard form of the differential equation for the confluent hy- 
pergeometric function. The solution with the required asymptotic behavior 
is [r(a)/r(c)]$(a, c; y), and we finally get 



8 Critical curves and bosonic fields (Coulomb gas) 

In the rest of this review I will provide a connection between SLE and a more 
traditional approach to critical 2D systems, namely, conformal field theory 
(CFT). In this Section we will se how critical curves can be described within 
a CFT of a scalar field. Closely related discussions have appeared before in 

Refs. IsniEZlEHl- 



8.1 Prom loop models to bosonic fields 

The relation between critical curves and operators of a boundary CFT is 
most transparent in their representation by a Gaussian boson field z) 
jl9l 1^2. 59, 60^. This representation is commonly known as the Coulomb 
gas method. Specifically, let us consider the 0(n) model on a honeycomb 
lattice. In the hope of describing the critical point by a local field theory, 
we need to have a description in terms of local weights on the lattice. 

To reproduce the partition function we randomly assign orientations 
to loops and then sum over all possible arrangements. The sum of weights 
for two orientations of every loop should give n. This is achieved by giving 
the local weight g^*^"'^/^ to each lattice site where an oriented loop makes 
right (left) turn. The weight of an oriented closed loop is the product of 
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all local site weights, and is equal to e^^'^°'" since for a closed loop the 
difference between the numbers of right and left turns is ±6. The sum over 
the orientations reproduces the correct weight n for an un-oriented loop if 
we choose 

n = 2 cos vreo- 

The range of —2 ^ n ^ 2 where the loops of 0(n) are critical can be covered 
once by eo G [0,1]- However, as we will see, to describe both the dilute 
and the dense phases we need to allow for a wider range cq S [—1, 1], with 
positive eo for the dense phase and negative eo for the dilute phase. 

For each configuration of oriented loops we can define a real height vari- 
able H that resides on the dual lattice and takes discrete values convention- 
ally chosen to be multiples of tt. To define H we start at some reference 
point where we set H = 0, and then every time we cross an oriented loop, 
we change H by zbvr depending on whether we cross the loop from its left to 
its right side or vice versa. Since the orientation of the loops was introduced 
randomly, the height function has to be compactified with radius TZ = 1: 

~ + 27r, (72) 

which means that the heights H and H + 2tt correspond to the same con- 
figuration of un-oriented loops. 

At criticality, the coarse-grained height function becomes a continu- 
ous scalar field (boson), believed to be described by the Gaussian action 
(g/^Tr) Jjy (Px (Viif)^, where the fluctuation strength parameter g is not yet 
determined. This can be done either by comparison with exact solutions 
of a related six-vertex model, or by an elegant argument due to Kondev 
and Henley |611 162j (which, unfortunately, only works in the dense phase). 
Here's the argument. 

If the system is defined on a domain with boundaries, some loops may not 
be counted with the correct statistical weight. For example, the difference 
between the numbers of left and right turns for a loop that wraps around a 
cylinder is rather than 6. Therefore, without modifications all such loops 
will be counted with a wrong weight 2 in the partition function . This is 
fixed by adding to the action a boundary term (ieo/27r) Jq^cU KH, where 
K is the geodesic curvature of the boundary. Each loop wrapped around 
the cylinder introduces an additional height difference AH = zLtt between 
the ends thus acquiring the correct weight. 

A similar situation occurs if the critical system lives on a surface with 
curvature, which microscopically can be viewed as existence of defects on 
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the honeycomb lattice (pentagons and heptagons correspond to positive and 
negative curvature, correspondingly). The correct weight for a loop that 
surrounds a region of non-zero curvature is obtained only if we include in 
the action the so-called background charge term (zeo/Svr) (fx RH, where 
R is the scalar curvature. 

Yet another necessary term in the action is the locking potential of the 
form A J (PxV{H) which would force the discrete values H = kir in the 
limit A oo. It must be, therefore, a 7r-periodic function of H, the most 
general form of it being V = X^^g^ ^.^q v^e^^^^ . Each term of y is a vertex 
operator whose dimension is j22j 

2 

Xk = -k{k - eo). 

g 

Most of these terms are irrelevant at the Gaussian fixed point, and we can 
ignore them. The most relevant term has /c = lifO<eo<l, and it has to 
be strictly marginal {xk = 2) in order to retain the conformal invariance of 
the action. This gives the relations 

eo = l — g, n = —2cosTTg. (73) 

In this case < 5 < 1, which is known to describe the dense phase of the 
0(n) model. In the dilute phase the second relation (|7!^ still holds [22]) 
but with 1 ^ g( ^ 2. This range is not possible to obtain from the previous 
argument since for — 1 < eg < we would need to pick k = —1 term as the 
most relevant, and it would still give us(7 = l + eo = l — |eo| <!• With some 
amount of hindsight we will assume both relations (|73)) to be valid for the 
whole range g G (0, 2] encompassing both the dense and the dilute phase. 
The point g = 1 separating the phases is somewhat special: there we need 
to keep both k = 1 and k = —1 terms in the locking potential since they 
have the same dimension. 

The failure of Kondev's argument in the dilute phase has a very signifi- 
cant geometric meaning. Namely, upon the coarse-graining the 0(n) loops 
become level lines of the bosonic field. However this identification can only 
be made for the dense phase, where the loops come close to themselves and 
each other on the lattice, translating to them becoming non-simple curves 
(with double points) in the continuum limit, resembling the traces of SLE 
with K > 4. The relation between critical lines and the bosonic field is quite 
different in the dilute phase, and this difference is related to quite a few sub- 
tleties in the treatment of both the dilute and the dense phases of a bounded 
system in the Coulomb gas formalism. For details see our paper |21j . 



53 



We now introduce the parametrization 



4 

g = -, 2 ^ K<oo, (74) 

K 

where k can be identified with the SLE parameter by comparing calculations 
of some quantity within the two approaches. A typical example is the dis- 
tribution of winding angles of critical curves on a cylinder, which is known 
through both the Coulomb gas method and SLE. Another good example 
is the multifractal exponents related to derivative expectations, which we 
compute in the Coulomb gas formalism in Section IHl Notice that k < 4 and 
K > 4 describe the dilute and the dense phases, correspondingly, while k = 4 
gives the point g = 1 separating the two phases. All this is quite consistent 
with the SLE phases determined in Section [6.21 

In the CFT literature it is customary to rescale the field (/? = \/2gH and 
make the coupling constant fixed 5ncw = 1/2, at the expense of varying the 
compactification radius of ip: 

n = vW^- (75) 

This is the normalization that we adopt from now on. For the rescaled field 
the action with all the terms becomes 

S[^] = — I (fx \{V^f + i2V2aQR:A + I dl Kip 

Jd 27r Jgij 

+ [ d'^xe'^''+^, (76) 
Jd 

where we use the notation 

2 \/k' 



2«o = 7=, a± = ao± \/ Oq + 1, 



a+ = ^/^/2, a- = -2/^. (77) 

Notice that oq is proportional to eo, and can be both positive and negative, 
its sign being different in the two phases of the loop model. 



8.2 Coulomb gas CFT in the bulk 

Consider now our bosonic theory on the infinite plane (the Riemann sphere) , 
dropping for now the boundary term in Eq. H76j) . The action S[(p] does not 
describe a free field beacuse of the presence of the locking potential. In 
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practice, however, this potential is always treated perturbatively, and any 
correlation function is expanded as 

oo 

{X)s = d^xi ... (i2x„(e*^"+^(^i) . . . e^^°+^(^")x), (78) 

n=0 ^' •' •' 

where (...) stands for correlators in the free theory with action 

So = ^ f cPx + i2V2aQR^] . (79) 

The neutrality condition discussed below makes sure that for a given oper- 
ator X, at most one term survives in the sum in Eq. (|78)) . The free action 
So is known to describe a CFT with the central charge 

= 1 - 2Aal = 1 - 3 ^ \ (80) 

which is the same as Eq. ()18() . The holomorphic part of the stress-energy 
tensor corresponding to the central charge (|80|) is 

T = -]^:{d^f:+iV2aod^V, (81) 

where d = d/dz, and semicolons stand for normal ordering. 

Notice that by the appropriate choice of metric, the curvature R may be 
made to vanish everywhere in the finite region of the plane. The curvature 
is then concentrated at infinity, and its effect is represented by insertion of 
a certain vertex operator (l^-2ao,-2oo ™ notation of Eq. (|51)) below) 
in the correlation functions, thereby changing the neutrality condition, see 
discussion below. This prescription, due to Dotsenko and Fateev PSl; allows 
to calculate correlators of primary fields (vertex operators) using simple free 
boson with c = 1 described by the action ((7^ but with ao = 0: 

So = ^ [ d^xiV^f. (82) 

In the complex coordinates z = x + iy,z = x — iy the field ip separates 
into the holomorphic and antiholomorphic parts: 

ip{z,z) = (l){z) + 4){z), 

and the basic correlators of these fields follow from ()82() : 

{^{z)4>{z')) = -log{z-z'), {mH^')) = -^og{z-z'), (0(z)0(z)) = 0. 

(83) 
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Primary fields in the theory (|79() are electric (vertex) and magnetic (vor- 
tex) operators, and their combinations also called vertex operators for sim- 
plicity (they all are implicitly assumed to be normal ordered): 

V,,,{z, z) = e*v^'^'^(^'^"), Vo,U^, z) = e-v^™^(^'^"), 

where we introduced the Cauchy-Riemann dual 

lf{z, z) = -i(t){z) + i(j){z) 

of the field f, as well as electric and magnetic charges e and m. A gen- 
eral vertex operator can also be written as a product of holomorphic and 
antiholomorphic components: 

Va,a{^^ ^) = Va{z)Va{z) = e^^'^'^i^') e'^'^^^^\ (84) 

where the holomorphic and antiholomorphic charges are: 

a = e + m, a = e — m. 

The holomorphic and antiholomorphic dimensions of the vertex opera- 
tors follow from the anomalous stress-energy tensor l|81j) : 

h{a) = a{a — 2ao) = h{e, m) = {e + m){e + m — 2ao), 
h{a) = a{a — 2ao) = h{e, m) = {e — m){e — m — 2ao). 

From this we see that a vertex operator is spinless (meaning that h = h) if 
either a = a or a = 2ao — Oi- In the first case the operator is purely electric 
(m = 0), and in the second case it can have an arbitrary magnetic charge, 
but the electric charge should be e = ao- We then introduce the notation 

V^''\z,z) = Va{z)V2ao-a{z). (85) 

Notice also a certain duality: the dimensions of the operators Va and 
V2ao-a are the same. This is consistent with the correlator 

{Va{z)V2a,^a{z')) = (,Z-z'r"'-. 

We see that the sum of the charges of the operators within the correlator is 
2ao, which is the negative of the background charge —2ao placed at infinity. 
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This is true in general: in the theory with a background charge correlators 
of vertex operators do not vanish only if the following neutrality condition 
is satisfied — the total sum of charges should equal to 2ao, in which case 
the chiral correlator is given by 

(^1)^42(^2) •..K.„(2n)) = n(^^-^i)'"'"'' E"* = 2«o. (86) 

i<j i 

Similarly, the correlator of vertex operators Va,cx{z,z) is 

(n^°-°»(^^'^*)> = ni^*-^il'°'°'' 5^ai = 2ao. (87) 

i i<j i 

While the global behavior of correlators of vertex operators is affected 
by the background charge, their local properties are completely encoded in 
the short-distance operator product expansions (OPE) 

VaAzi)Va,{z2) = (21 - Z2)'"^"'Ki+a2(^2) + ■ • • , (88) 

where the dots stand for subleading terms. 

Finally, let us mention that in CFT literature it is customary to label 
holomorphic charges and weights by two numbers r, s according to 

ar,s = ^(1 - r)a+ + ^(1 - s)a_, 

K 2 \2 2 \2- 

2~ ~ 7^ 



arAar,s-2ao) = -[[r—-s-^) "( 
{rK - 4s)^ - (n- 4)^ 



16k 



(89) 



We stress that we use it as a shorthand notation and do not impose any 
restrictions on r or s. The holomorphic primary field of weight /i^^s is denoted 
by ^pr,s{z). In the notation (jHH, this is the field Va^ ,.{z). The corresponding 
spinless bulk operator V^"'''''\z, z) will be denoted as ^lJr,s{z,z). 



8.3 Coulomb gas CFT in the upper half plane 

We now consider modifications to the Coulomb gas description that result 
when the boson lives in a bounded region. Due to conformal invariance we 
may choose the simplest possible case: the upper half plane H, and most 
formulas will be written for this case. For basics on boundary CFT see Refs. 

[iniiSlEll- 
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First we note that since the boson field is defined by the orientation 
of the loops, it is a pseudoscalar, meaning that it changes sign under parity 
transformation, or reflection in a boundary [z ^ z for the upper half plane). 
This implies that we must impose the Dirichlet boundary condition 

dMoD = 0, (90) 

where the derivative is taken along the boundary (for EI the boundary is the 
real axis z = z). Each boundary is then a level line of in correspondence 
with lattice loops. 

The Dirichlet boundary condition glues together the holomorphic and 
antiholomorphic sectors of the theory in a way that is easiest to describe 
in terms of the image charges, see, for example. Chapter 11 in Ref. jl9j . 
The dependence of correlators of primary flelds on the antiholomorphic co- 
ordinates Zi in the upper half plane can be regarded as the dependence on 
holomorphic coordinates in the lower half plane, after the parity trans- 
formation is performed. In our case the parity transformation for the chiral 
boson is simply 



which gives the following prescription for bulk vertex operators in the upper 
half plane: 

y^^(^Z,z) = e*V^^'^(^'^) = (,iV2e[<j,{z)+4>{z)] _^ ^iV2e<t>(z)^~iV2e<t>(z*)^ 
VQm{z,z) = e"^™*^*^^'^) = e*^'"['^(^)~'^(^)] ^iV2m<j)(z) ^i\/2m<l){z*) ^ 
Vem{z,z) giV2(m+e)(l){z) ^iV2{m~e)(j)(z*) ^ 

Vaa{z,z) = f^y^'^<t>{^) ^i-y^^'i'i^) giV2a(j){z) ^-iV2acj>{z*) ^ ^g-^^ 

The right hand sides of these equations should be viewed as products of two 
holomorphic operators. 

This situation can be summarized by saying that under reflection the 
electric charges change sign, while the magnetic ones do not. This implies, 
in particular, that on the boundary (the real axis) only magnetic operators 
survive, electric ones being rendered trivial by the Dirichlet boundary con- 
dition. Indeed, as z and z* both approach a point x on the real axis, the 
bulk operator Ve^m{z-,z) (with any electric charge e) reduces to 

y (2m) ^^-j _ gi2y2m(/)(a:) ^ (92) 
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Such boundary operator is characterized only by one weight which is the 
same as the holomorphic weight h{0, 2m) of a bulk operator. This situation 
can also be described by saying that the fusion on the boundary of a holo- 
morphic operator with its image produces a magnetic boundary operator. 



8.4 Creation of critical curves 

As we have seen, in a microscopic description a critical curve starting at 
a boundary is created by a change in boundary conditions. In the effec- 
tive language of CFT such change is implemented by insertion of a certain 
operator at a point on the boundary. 

Microscopic definition of height function H implies that if we have n 
critical curve with the same orientation that start on a boundary at some 
point, the boundary values of the field ip on the two sides of the point should 
differ by zbnvrT^. The curves should be oriented in the same way to prevent 
them from reconnecting with each other. For a system on the UHP, n such 
curves starting at the origin are then created by the following boundary 
condition: 

, . / n7^7^, X ^ 0, , . 

= 1 0, X > 0. (9^) 

The operator creating such a jump in the value of at x = is a cer- 
tain magnetic boundary operator. To find out what it is, we first consider 
magnetic operators in the bulk. 

Note that the insertion of a magnetic operator in the bulk creates a 
vortex in the field (p. Indeed, the bulk OPE l|88|) gives 

(2; z' \ 2e™ _ 
Ve+m{z')Ve-m{z') + ■ . . 
Z — Z J 

= e^^^'"^'-s(^--')Fe+m(^')V'e-™(^') + • • • • 

This means that when z goes around z', the field p changes by 4\/2vrm, and 
hence, a discontinuity line with this jump arises, see Fig. 1131 If this vortex 
corresponds to a star of n critical curves joined at the point z', the change 
in (p should be equal to mrTZ, and then the discontinuity is not physical due 
to the compactification of (p. This gives the magnetic charge of the bulk 
curve creating operator: 

\/2 ^ n n 
m = — nTc = — = = ——a-, (94 
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Figure 13: A bulk magnetic operator Vo,m creates a vortex configuration of 
the field (p. The field changes by 2\/87rm when going around Vo,m- 

where we used the value (|75|) of the compactification radius and the def- 
inition H77|) of a_. In order to be spinless (otherwise the operator would 
transform under rotations, giving non-trivial dependence on the winding 
number of curves), the bulk curve creating operator should also have elec- 
tric charge qq. 

We have, therefore, found the holomorphic charge of the bulk curve 
creation operator to be 

n 

a = "0 - = ao,n/2 

in the notation of Eq. (|89j) . The operator itself is then V'o,n/2('2, -z), and its 
holomorphic weight is 

_ 4n2 - (k - 4)2 

"'0,n/2 — 77 • 

In particular, a single critical curve going through a point z is created by the 
operator •i/'o,i(-2, ^) with the holomorphic weight /io,i = (8 — k)/16. Notice 
that this weight is related to the fractal dimension of the critical curve by 
df = 2 — 2/io,i, see Ref. IHSl for details. 

Now we go back to the boundary. According to Eq. H91() . as z ap- 
proaches a point X on the real axis, the bulk operator Ve,miz,z) (with any 
electric charge e) reduces to e*2v^"*<^(^) . In this process the two sides of 
the discontinuity line become parts of the real axis separated by x (see Fig. 
I14() . Thus, when we go from one side of x to the other along a semicircle, 
the field changes by the same amount as in making a full circle around a 
bulk magnetic operator. Then to create the boundary condition (|93|1 . and, 
correspondingly, n critical curves starting at the origin, we need to insert 
there the magnetic operator 
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V, 



0,m/2 



V8 







Figure 14: A boundary magnetic operator is obtained as the boundary limit 
of the bulk magnetic operator V(j „^/2- 



with the magnetic charge determined by the condition \/87rm = mrTZ, which 
gives 



I n n 

m = —=niZ = —= = a_. 

\/k 2 



(95) 



In the notation of Eq. (|S^ this is ai^n+i-, so the boundary curve creating 
operator is ipi^n+i{x) with dimension 



h 



l,n+l 



2n^ + n(4 - k) 
2k 



In particular, a single curve is created by the insertion of V'i,2(3;) with the 
dimension 

6 — K 



hi: 



2k 



(96) 



which appeared before in the discussion of restriction property of SLE, see 
Eq. ((SI. 



9 Harmonic measure of critical curves 

Harmonic measure is a useful quantity describing geometry of complicated 
plane domains [661 16 7j . In the following sections I define it and the related 
spectrum of multifractal exponents, and then show how to compute these 
exponents for harmonic measure on critical curves using CFT. 
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9.1 Definitions of hcirmonic measure 

Let D he a domain (an open connected subset of the complex plane), dD 
its boundary, and z e D. 

Harmonic measure in D from z, denoted as uJij{z,T) where F G dD, is 
a probability measure on dD, which can be defined as the probability that 
the standard two-dimensional Brownian motion Bt that starts at z hits dD 
in a given portion T C dD of the boundary: 

Here td is the escape time from D, that is, the first time when the Brownian 
motion Bt hits the boundary dD. 

Harmonic measure u)d{z, F) can also be characterized as the unique har- 
monic function (solution of the Laplace's equation) u{z) in D with the 
Dirichlet boundary conditions 

«(C) = i, CeF, u{0 = o, c^r. 

Here wc will only be interested in harmonic measure from infinity of 
the domain D exterior to a closed curve 7. In this case we will denote it 
simply by ci;(F). Harmonic measure ci;(F) has an electrostatic interpretation. 
Imagine that D is a charged metallic cluster with the total charge one. The 
charge is concentrated on the boundary dD. Then a;(F) is the charge located 
on the portion F of OD. 

Harmonic measure is conformally-invariant: if Z)' is a conformal 

map that is continuous and one-to-one on D U dD, then 

u;D{z,T)=u;D'{f{z),f{r)). 

9.2 Moments of heirmonic measure and multifractal spec- 
trum 

Consider a closed curve 7. One can cover the curve 7 with discs B{zi,r) of 
radius r centered at some points -^i G 7 {zi form a discrete subset of 7). Let 

p{zi,r) = u{'jn B{zi,r)) 

be the harmonic measure (from infinity) of the portion of the curve covered 
by the disc B{zi, r). Then we consider the moments 

N 

Mh = ^pizi,rf, (97) 

i=l 
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Figure 15: A curve covered by discs. 



where /i is a real power, and N is the number of discs needed to cover 7. 
As the radius r gets smaller and the number of discs gets larger, these 
moments scale as 



The size L (diameter) of the curve is used to make the right hand side of 
this equation dimensionless. 

The function t(/i) is called the multifractal spectrum of the curve 
7. This function encodes a lot of information of the curve 7. It also has 
some simple properties. First of all, since all < p(zi,r) ^ 1, the moments 
M/j are well defined for any real h, and the function t(/i) is non-decreasing: 
r(/i) ^ T{h') for any h < h' . Secondly, if /i = 1, the sum in is equal to 
the total charge of the cluster, and therefore does not scale with r, producing 
the normalization condition 



Third, if we set h = 0, Mq is simply the number A^ of discs of radius r 
necessary to cover the curve 7, so that by definition the fractal (Hausdorff) 
dimension of 7 is 



If the curve 7 were smooth, we would have a simple relation t(/i) = h — 1. 
For a fractal curve one defines the anomalous exponents 6{h) by 




(98) 



r(l) = 0. 



r(0). 



T{h) =h-l + 6{h). 
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Also, the generalized multifractal dimensions of a fractal curve 7 are defined 
as D{h) = T{h)/{h — 1) (so that D(0) = df). A non-trivial theorem due to 
N. Makarov [HHl EH states that 



Dil) 



r'(l) = 1. 



If the curve 7 is a member of an ensemble of curves generated in some 
random way, the moments become random variables, and we can study 
their distribution functions. If the distribution functions are narrow, then 
the mean moments M/j, where the over line denotes the ensemble averaging, 
characterize them well. This is usually the case for \h\ that is not too large. 
In this case the mean and the typical moments scale in the same way, which 
is a statement about self- averaging of the moments. This is what we assume 
here. 

In this situation it is also natural to assume some sort of ergodicity, 
meaning that the summation over the points Zi in Eq. (|97() for a typical 
curve is equivalent to the ensemble average. Hence we can write 



where now the harmonic measure p{zQ,r) is evaluated at any point zq G 7. 
We define the local multifractal exponent T{h) at a point zq by 



Similar to t(/i), for a smooth curve we have r(/i) = h, so in general we define 

(2) 

the local anomalous exponents A^^[j,(/i) by 



The reason for the superscript (2) and the subscript "bulk" will become 
clear in the next subsection. 

^ ('2) 

It is obvious from the definitions that r(0) = A^^[j,(0) = 0, and we 
deduce simple ergodicity relations 





(99) 



?(/!)= /i + A(l(/i). 



T{h) 

df 



t(/i)-t(0), 

m-m, 




r{h) 

m 



?(/i)-f(l), 

Af,l(/^) - a[,1(1) 



(100) 
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9.3 Critical curves and uniformizing maps 

So far we have considered arbitrary closed curves. An example of such a 
curve is the exterior perimeter 7 of a critical cluster. One can imagine that 
the cluster is made of a conducting material and carries the total unit electric 
charge. The harmonic measure of any part of 7 is then equal to the electric 
charge of this part. Since exterior perimeters are always dilute curves in the 
sense of Section IH?T] f see also discussion of duality around Eq. ^^), in the 
remaining part of this paper we will always assume k ^ 4. 

The critical clusters and their boundaries appear as members of statisti- 
cal ensembles, which is the situation suitable for local multifractal analysis 
of the previous section. We then will pick a point of interest zq on the curve 
7 and consider a disc of a small radius r <^ L centered at zq. It surrounds 
a small part of 7, and we will study the mean moments of the harmonic 
measure p{zo,r) and their scaling as in Eq. 

There are a few generalizations of the simple closed critical curve con- 
sidered above. First of all, the curve 7 need not be closed or stay away from 
system boundaries. If 7 touches a boundary, we can supplement it with 
the image 7 (reflected in the boundary) and take the union 7 U 7 to be the 
charged conducting object. The electrostatic definition of p{zQ,r) can be 
naturally extended to the cases when zq is the endpoint of n critical curves 
on the boundary or in the bulk. If n is even, the latter case can be also seen 
as n/2 critical curves passing through zq. In particular, n = 2 corresponds 
to the situation of a single curve in the bulk, considered above. 

When Zq is the endpoint of n critical curves on the boundary or the bulk 
we define the corresponding scaling exponents similar to Eq. (|99() : 

p{zo,r)h ~ r'^+Af-'W^ p{zo,r)'' ~ r'^+^iuLW. (loi) 

In the case of a single curve we will drop the superscript, so, for example, 
A{h) = A(i)(/i) is the same exponent as obtained in Eq. (|71|) . 

These exponents were first obtained by means of quantum gravity in 
[21 El • For a critical system with parameter k the results read 



^ K - 4 + v/(k - 4)2 + 16K/i ^ Vi-c + 24/t- 

A(")(/i) = nA(/i), (103) 
a(1(/^) = --:^ + (^ + ^) (- - 4 + V(-- 4)2 + 16k/.). (104) 



2 V16 4k 

Remarkably, A(/i) is the gravitationally dressed dimension h, as given by 
the KPZ formula of 2D quantum gravity [^jZOl ■ Starting in the next section. 
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Figure 16: The uniformizing confer mal maps for various cases considered. 
The dots denote that points where the electric field is measured. 

we will show how to obtain these exponents in the framework of Coulomb 
gas CFT, where they are also written in a more transparent way. 

The basic property that allows ro calculate the multifractal exponents 
using CFT is the conformal invariance of the harmonic measure. Let us 
consider a conformal map w{z) of the exterior of 7 to a standard domain. 
Usually we will choose the upper half plane but sometimes the exterior of 
a unit circle is more convenient. We normalize the map so that the point 
of interest zq is mapped onto itself, choose it to be the origin zq = and 
demand that at infinity w{z) = z + o{z). Examples of w{z) for several cases 
are shown in Fig. 1161 

The scaling of it;(z) near the origin is directly related to that of the 
harmonic measure. Indeed, since p(0, r) is the charge inside the disc of radius 
r, by Gauss theorem it is equal to the flux of the electric field through the 
boundary of this disc. This, in turn, should scale as the circumference of 
this disc times a typical absolute value of the electric field at the distance r 
from the origin, i.e. |^i;'(r)|. This leads to scaling relation 



The relation of the scaling of the harmonic measure and the derivative 
of a uniformizing map allows for further generalizations. Namely, we can 



p(0, r) ~ r|i(;'(r)| 



which allows to rewrite the definitions ()1U1() as 





(105) 
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measure the electric field in more than one point. Close to the origin n 
curves divide the plain into n sectors in the bulk and n + 1 on the boundary. 
Then we can study objects like 

. . . \w' {zn+i)\^'^+^ (boundary), 
\w'{zi)\f'^ ...\w'{zn)\^" (bulk), (106) 

where no two Zj's lie in the same sector. The case when the electric field is 
not measured in some sectors is done by setting hi = in them. We will see 
how to express these quantities as CFT correlation functions. In the case 
when Zi are all at the distance r from the origin {zi = re*^% 6i = const), 
these averages scale as 'r^''"'Hhi,-hn+i) g^j^j j,A^^j^(/ii,...ft,„) ^j^^j^ ^j^g higher 
multifractal exponents (HI 

n+l n+l 

^A(")(/i,) + ^j;A(/^.)A(/i,), (107) 

i=l i<j 
n n 

E^bulk(^^) + -^Y.H^^)Hh,)■ (108) 

i=l i<j 



A^^\hl,...hn+l) = 

Airi(/ii,...M = 



9.4 Derivative expectations and CFT in fluctuating geome- 
try 

Here we begin the Coulomb gas derivation of results ()1U2I - [TU11 11071 11U8() . 

It is easiest to start with a point where a single curve 7 connects with the 
system boundary. We assume that the critical system occupies the upper 
half plane, so that the real axis is the boundary. 

The partition function Z{0,L), restricted to configurations that contain 
a curve 7 connecting the points and L, is given by the correlator of two 
boundary curve creating operators, in this case ipi^2 (see section (|8.4|) ): 

= (V'l,2(0)^l,2(i)>H, 

where Z is the unrestricted partition function. This correlation function can 
be computed in two steps. In the first step we pick a particular realization 
of the curve 7. Within each realization, it is the boundary separating two 
independent systems — the interior and the exterior of 7. In both these sys- 
tems we can sum over microscopic degrees of freedom to obtain the partition 
functions Z™* and respectively. These are stochastic objects that de- 
pend on the fluctuating geometry of 7. In the second step we average over 
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the ensemble of curves of 7. We thus obtain 



Z{0,L) = ZfZ^^K 

Next, we insert an additional boundary primary operator Oh{r) of di- 
mension h close to 0, and another one O/i (00) at infinity. The first one serves 
as a "probe" of harmonic measure, and the second is necessary to ensure 
the charge neutrality. We thus consider the correlation function 

(V'i,2(0)O/,(r)Vi,2(^)Oh(oo))H. (109) 

Since we are only interested in the r-dependence of the correlation function, 
we can fuse together the distant primary fields: ipi^2iL)Oh{oo) ^'(oo). 
We therefore consider the r-dependence of a 3-point function 

(^i,2(0)O/,(r)^'(oo))H, (110) 

and show that it yields the statistics of the harmonic measure. 

Decomposing the upper half plane into the exterior and the interior of 7 
as before, we can rewrite ()1U9() as the average over the fluctuating geometry 
of 7: 

(0;,(r)0^oo))-tZij'tZ-t. (Ill) 

Here the domain of the definition of the correlation function of primary 
fields is the exterior of 7. This correlation function is statistically inde- 
pendent from the other two factors in the numerator of (|11H) in the limit 
r <C and we are left with the correlation function {Oh{r)Oh{oo))'^^ of 
two primary fields of boundary CFT, further averaged over all configurations 
of the boundary 7. This average is proportional to the 3-point correlation 
function ^TU^. 

Now we apply the conformal transformation w{z) which maps the exte- 
rior of 7 onto the upper half plane. Being a primary operator of weight h, 
Oh{r) transforms as Oh —>■ \w'{r)\^Oh{w{r)), while Oh{oo) does not change 
because of the normalization of w{z) at infinity. The transformation relates 
the correlation function in the exterior of 7 to a correlation function in the 
upper half plane: 

{Oh{r)Oh{oo))f = \w\rt{Ohiwir))Oh{oo))M. (112) 
The latter does not depend on r and can be neglected. 
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Summing up, we obtain a scaling relation between the moments of the 
harmonic measure near the boundary and correlation functions of primary 
boundary fields [HS] : 

K(OF « (O/^WV'l,2(0)^(oo))H, r < |L|. (113) 

The primary field ^'(oo) should be chosen in such a way as to render the 
correlation function non-zero. The choice is made unique by picking the 
conformal block which satisfies simple physical condition A(0) = 0. The 
r-dependence of the correlation function (|113() is found from the OPE of the 
fields Oh{r) and ^1,2(0): 

00 

Oh Wi,2 (0) = ^r^'^- (0). (114) 

k=l 

The exponent A(/i) is then identified as the lowest power A^. such that 
($(fc)(0)^(cx))) / 0. 

Several remarks are in order. As presented, this argument produces the 
scaling exponent A(/i) for a single curve on the boundary. But it can be 
easily modified for studying the scaling behavior in all other cases. The case 
of n curves starting from a point on the boundary is obtained by simple 
replacement of the curve creating operators: V'1,2 — *■ V'i,n+i (see Section 

Also, the argument can be repeated for the case of 7 connected to the real 
axis only at one point (as in SLE). In this case no separation in two systems 
is necessary. Finally, nothing compels us measure the electric field on the 
real axis. We could take instead a bulk primary field 0/i/^/j/(z, z), where 
\z\ = r. The weight h' should be chosen such that when the holomorphic 
part Oh'{z) is fused with its image Oh'{z*), the boundary field Oh{ ^\^ ) 
with dimension h is obtained, similar to Eq. (|92)) . This fusion will be used 
below. 



9.5 Calculation of boundary multifractal exponents 

In practice we view both Oh and Vi,2 in Eq. (|114|) as boundary vertex 
operators and l/('^i>2)^ gee Eq. ()92() . with charges 

a_ 1 



Oh = ao - \/ + h, ai^2 /-- 
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The leading term in the OPE of these two operators corresponds to simple 
addition of charges, see Eq. (|88() . Hence, the scaling relation ()113|) immedi- 
ately gives the result H1U2() written in a compact and suggestive form: 

A(^) = 2ai,2ah. (115) 

It is interesting that written in this form, the KPZ formula for gravitationally 
dressed dimensions amounts to OPE of vertex operators in a simple Coulomb 
gas OFT, without any quantum gravity. 

An immediate generalization to the statistics of harmonic measure of 
n curves reaching the system boundary at the same point is obtained by 
replacing ^1^2 V'i,n+i- Since ai^n+i = —na-/2 = nai,2, this immediately 
leads to 

aW(/i) =2Qi,„+ia;, = nA(/i), 

which is the same as Eq. ()1U3() . 

To calculate the higher boundary multifractal exponents, we consider n 
non-intersecting critical curves growing from the origin on the boundary (the 
real axis). It will be convenient to assume that the curves end somewhere 
in the bulk thus forming a boundary star (e.g. the third picture in Fig. 116(1 . 
Let w{z) be the conformal map of the exterior of the star to the upper half 
plane with the usual normalization w{z) = z + o{z) at z ^ 00. 

We want to find the scaling of the average 

\w'{zi)\h^ ...\w'{Zn+l)\^"+\ 

where zi are all close to the origin, no two of them lying in the same sector. 
The latter condition will be automatically satisfied in the subsequent calcu- 
lation due to the following: if in a particular realization two points Zi and 
Zj happen to be in the same sector, then w{zi) — w{zj) — > as — Zj — > 0, 
but if they lie in different sectors, w{zi) — w{zj) remains large in the same 
limit. 

Since n curves starting from the origin on the boundary are produced 
by the operator V'i,n+i(0), we now consider a boundary CFT correlation 
function with several "probes" of the harmonic measure: 

n+l 

C = {l{O^,,^,,:^{zi,z,)i;i,n+i{0)-^'{oc))^. (116) 
1=1 

The primary field at infinity represents the fusion of all fields far from the 
origin and should be chosen by the charge neutrality condition. As before. 
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this correlation function is equal to the statistical average of a certain corre- 
lator in the fluctuating domain, and we further apply the uniformizing map 
w(z) to transform this domain into the UHP: 



i i 

Unlike Eq. ()112() , the correlator under the average cannot be neglected since 
it does depend on the short scale r, as we shall see soon. 

The correlator C can now be evaluated in two ways. As before, we view 
the primaries Oj^i i^i{zi^Zi) as vertex operators with charges 



a'i = UQ- ^al + h[. (118) 

Then, using the prescription (|9H) . we can rewrite C in Eq. (|116j) as a full 
plane chiral correlator, which then is evaluated using Eq. H86() : 

C = ( n O/.: {Z,)0h'^ «) Vl,n+1 (0)^(00)> 
i 

i i<j i,j 

When all Zj are at the same distance r from the origin, the last expression 
scales as 

On the other hand, we can evaluate the correlator that appears inside 
the average in Eq. (|117|) in the same way: 

i i 

OC X{{w{zi) - w*{zi)f<" W \w{zi) - «;(z,)l'"^"^- H (^(^^) " ^(4))'"^"^- 

i i<j 

We specifically separated the diagonal (i = j) terms, since only they con- 
tribute to the necessary short-distance behavior. All the other term insure 
that the realizations of the curves in which any two points Zi end up in 
the same sector are suppressed (since the distances w{zi) — w{zj) are then 
small), and we can consider only the case when all w{zi) are far apart. Then 
the relevant short-distance dependence of Eq. (|117|) is 



C oc J] \w'{zi)\''^^{w{zi) - w*{z,)y<' 
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Insofar as the scaling with r is concerned, we further approximate w{zi) — 
w*{zi) ~ |zj| [^'(zi)! ~ rlw'^Zi^l. This gives 

C oc r^^^''''l[\w>{zi)\^f^^+^<\ (120) 

i 

The exponents in the last factor 

2/i'j + 2a'j^ = 2a^(a^ — 2ao) + lol-^ = 2a'j(2a^ — 2ao) = /io^ = hi 
are the dimensions of the boundary operators with charges 



appearing in the OPE of two chiral operators with charges a[. 
Finally, comparing Eqs. pi9j) and (|12()j) . we get the result 

\w'{zi)\''^ ...\w'{Zn+l)\^-+^ (X r^'-K/^i.-./^n+i)^ 

with the higher multifractal exponent 

n+l n+1 
i=l i<j 

which is the formula 1)107^ . 

9.6 Calculation of bulk multifractal exponents 

Calculation of bulk multifractal behavior is done in much the same way as on 
the boundary, so we go straight to the general case of higher bulk exponents. 

Let the critical system, occupying the whole complex plane, be restricted 
to having n critical curves growing from a single point, in which we place 
the origin z = 0. We will assume that z = is the only common point 
of these curves, since the local results around this points are unaffected by 
the curves' behavior at large distances. We define the conformal map w{z) 
of the exterior of the "star" to the exterior of a unit circle with the the 
normalization w{z) = z + o{z) at z — > oo. 

Close to the origin the curves divide the plane into n sectors. We consider 
a quantity 

\w'{zi)\^^ . . . \w'{Zn)\^'^, 
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where Zi are points close to the origin, no two of them lying in one sector. 
As before, if two points Zi and Zj happen to be in the same sector, w{zi) — 
w{zj) — > when Zi — Zj 0, but if they lie in different sectors, wlzi) — w{zj) 
remains large. 

Since n curves starting from the origin in the bulk are produced by the 
operator ^^0,11/2(0)) we introduce a CFT correlation function 

n 

Cbulk = {Y[Oh>^^{zi,Zi)%^^,2{Q)^ {00)), (121) 

i=l 

where, as before, h'- is the weight of a primary field such that the result of 
its fusion with its image has the weight hi: 

ai = 2a'^. 

Proceeding exactly as in the boundary case, we first rewrite the correlator 
Cbuik as the ensemble average of another correlator in the exterior of the 
star of critical curves. Then we map that exterior to the exterior of a unit 
disk C \ D (see the second picture on Fig. [T^ : 

Cbuik=ni^'(^oi'''Kn^''^'^^(^(^*)'^(^*))^(°°)>c\D- (122) 

i i 

Next, we evaluate Cbuik as defined in Eq. H121|) . using Eq. (|S7|) : 

Cbuik = n n - ^ii^"^"' °^ r^"°'"/^^«"^+^^»<^"^"^'. 

i i<j 

Alternative evaluation starting from Eq. (|122|) gives the same result as Eq. 
(|12UI) . Combining the two results for Cbuik) we obtain 

with the higher bulk exponent 

n n 
i=l i<j 

where 

^buik('^) = 2ao,n/2ah - 2«/i = (2ao,n/2 - aoja/i - 

is the scaling exponent of a single |w;'(2;)|'^ in the presence of n critical curves 
in the bulk. These are the results quoted in Eqs. (|1()4[ ll()<^ . 
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10 Omitted topics: guide to the literature 



The current literature on SLE and related subject is already quite large. In 
this paper I had to omit many interesting and important topics. In this 
section I simply list the topics and give appropriate references. 

For the overall logic of this paper the biggest omission is the discussion of 
the relation of SLE and GET through the identification of GET correlators 
and SLE martingales. This identification was established and developed by 
Bauer and Bernard (see review and references there). We have further 
developed this correspondence |23 , showing that one can recover all familiar 
objects of the Goulomb gas GET, such as bosonic field, its current, vertex 
operators, and the stress-energy tensor, by focusing on SLE martingales. 
Alternative and independent versions of SLE-GFT correspondence and gen- 
eralizations were given by Eriedrich, Kalkkinen and Kontsevich, see Refs. 

Ghordal SLE considered in this paper has been generalized in many 
ways. First of all, one can define SLE in other simply-connected geometries 
than that of the UHP. The corresponding processes are known as radial 
SLE [II1IZ2], whole plane SLE [H], dipolar SLE [HSIIZSI. Ah these variants 
happen to be closely related [71| . 

Secondly, one can consider SLE in multiply-connected domains including 
arbitrary Riemann surfaces |^ EH ESI EOl HB HSl HHj ) though there is some 
amount of arbitrariness involved in the definition, since in this setting the 
conformal type of the domain changes during the evolution (one moves in 
the moduli space) pUl. 

Third, there is a way to modify the dynamics of the growth of a random 
curve by including certain moving points ("spectators") on the boundary 
of the domain, that influence the evolution by supplying a drift term in 
the (analog of) Loewner equation. This generalization is known as SLE^^p, 
where p stands for a vector of parameters describing the coupling of the 

"spectators" gHlEniinZlESIHIllHlllHSllHlI- 

The fourth generalization allows for multiple curves to grow simultane- 
ously. This is called multiple SLE |HS1 Ei EZl EHl ES EOl • 

We can also generalize SLE by dropping the demand that the forcing 
stochastic function be continuous, but keeping the requirement of station- 
ary and statistically-independent increments. This leads to a much broader 
class of forcing processes including, in particular, the so called Levy pro- 
cesses. This generalization might be a useful description of tree-like stochas- 
tic growth piPIT] . 

Still another generalization is to combine the evolution of conformal 
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maps with some stochastic process in a Lie algebra or some other algebraic 
structure. This leads to generalized SLE processes describing CFT with 
additional symmetries, such as Wess-Zumino models [nUESl- 

Recently, a few papers appeared that used SLE as a tool to probe con- 
formal invariance in systems that are not described by traditional statistical 
mechanics models. A remarkable example is Ref. |94| which numerically 
demonstrated that zero vorticity lines in highly developed 2D turbulence 
are SLEg with high accuracy. Similar conclusions were presented for do- 
main walls in spin glasses |95j and nodal domains of some chaotic maps 

m- 

Finally, I would like to mention that there is generalization of Loewner 
equation that describes evolving 2D domains which may grow with a speci- 
fied rate at every point on the boundary. These are called Loewner chains 
(see Refs. jSl ^ for a review of this enormous field in relation to SLE) and 
can describe various non conformally-invariant growth processes such as 
Laplacian growth, diffusion-limited aggregation, dielectric breakdown, etc. 
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